CONVERGENCE RATES FOR DISPERSIVE APPROXIMATION 
SCHEMES TO NONLINEAR SCHRODINGER EQUATIONS 
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Abstract. This article is devoted to the analysis of the convergence rates of several nu- 
merical approximation schemes for linear and nonlinear Schrodinger equations on the real 
line. Recently, the authors have introduced viscous and two-grid numerical approximation 
schemes that mimic at the discrete level the so-called Strichartz dispersive estimates of the 
continuous Schrodinger equation. This allows to guarantee the convergence of numerical ap- 
proximations for initial data in L^(R), a fact that can not be proved in the nonlinear setting 
for standard conservative schemes unless more regularity of the initial data is assumed. In 
the present article we obtain explicit convergence rates and prove that dispersive schemes 
fulfilling the Strichartz estimates are better behaved for //"(R) data if < s < 1/2. In- 
deed, while dispersive schemes ensure a polynomial convergence rate, non-dispersive ones 
only yield logarithmic decay rates. 



(1.1 
and 
(1.2 



1. Introduction 

Let us consider the linear (LSE) and the nonhnear (NSE) Schrodinger equations: 

iut + dlu = 0, X e M, t / 0, 
n(0, x) = (p{x), X € M 



iut + d^u = f{u), X G M, t / 0, 
u{0,x) = ip{x), X G R, 

respectively. 

The hnear equation (jl.ip is solved by u{x, t) = S{t)(f, where S{t) = e'*^ is the free 
Schrodinger operator and has two important properties. First, the conservation of the L^- 
norm 

(1-3) lk(*)llL2(R) = I|¥'IIl2(k) 

which shows that it is in fact a group of isometrics in L^(M), and a dispersive estimate of the 
form: 

(1-4) \SitMx)\ = |n(t,x)| < X G E, t / 0. 

The space-time estimate 
(1-5) II'S'(-)V5||l6(r,l6(k)) < C||v?||2,2(i.), 

due to Strichartz f29], guarantees that the solutions decay as t becomes large and that they 
gain some spatial integrability. 

Inequality (jl.Sp was generalized by Ginibre and Velo [llj. They proved: 

(1-6) \\S{-)y^\\Li{R,L-{R)) < Ciq)\\ip\\L2(K) 
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for the so-called 1/2- admissible pairs {q,r). We recall that the exponent pair {q,r) is a- 
admissible (cf. |23j) ii 2 < q,r < oo, {q,r,a) ^ {2,oo, 1) and 



We see that (jl.Sp is a particular instance of ()1.6p in which a = 1/2 and g = r = 6. 

The extension of these estimates to the inhomogeneous linear Schrodinger equation is due 
to Yajima ^32j and Cazenave and Weissler [6j. These estimates can also be extended to a 
larger class of equations for which the Laplacian is replaced by any self-adjoint operator such 
that the L°°-norm of the fundamental solution behaves like [23]. 

The Strichartz estimates play an important role in the proof of the well-posedness of the 
nonlinear Schrodinger equation. Typically they are used for nonlinearities for which the 
energy methods fail to provide well-posedness results. In this way, Tsutsumi [31] proved the 
existence and uniqueness for L^(R)-initial data for power-like nonlinearities F{u) = \u\Pu, in 
the range of exponents < p < 4. More precisely it was proved that the NSE is globally well 
posed in L°°(M, L^(M)) n L^^^(M, L^(W)), where {q,r) is a 1/2- admissible pair depending on 
the exponent p. This result was complemented by Cazenave and Weissler [7] who proved the 
local existence in the critical case p = 4. The case of iJ^-solutions was analyzed by Baillon, 
Cazenave and Figueira [T], Lin and Strauss [23], Ginibre and Velo [9i|10], Cazenave [1], and, 
in a more general context, by Kato [2H I22j. 

This analysis has been extended to semi-discrete numerical schemes for Schrodinger equa- 
tions by Ignat and Zuazua in [T7j, [18], |20j . In these articles it was first pointed out that 
conservative numerical schemes often fail to be dispersive, in the sense that numerical solu- 
tions do not fulfill the integrability properties above. This is due to the pathological behavior 
of high frequency spurious numerical solutions. Then several numerical schemes were devel- 
oped fulfilling the dispersive properties, uniformly in the mesh-parameter. In the sequel these 
schemes will be referred to as being dispersive. As proved in those articles these schemes may 
be used in the nonlinear context to prove convergence towards the solutions of the NSE, for the 
range of exponents p and the functional setting above. The analysis of fully discrete schemes 
was later developed in [TJ] where necessary and sufficient conditions were given guaranteeing 
that the dispersive properties of the continuous model are maintained uniformly with respect 
to the mesh-size parameters at the discrete level. The present paper is devoted to further 
analyze the convergence of these numerical schemes, the main goal being the obtention of 
convergence rates. 

Despite of the fact that non-dispersive schemes (in the sense that they do not satisfy the 
discrete analogue of (jl.Sp ) can not be applied directly in the L^-setting for nonlinear equations 
one could still use them by first approximating the L^-initial data by smooth ones. This paper 
is devoted to prove that, even if this is done, dispersive schemes are better behaved than the 
non-dispersive ones in what concerns the order of convergence for rough initial data. 

The main results of the paper are as follows. In Theorem 13.11 we prove that the er- 
ror committed when the LSE is approximated by a dispersive numerical scheme in the 
L''(0,T; /'"(/iZ))-norms is of the same order as the one classical consistency-|-stability analysis 
yields. Using the ideas of [3], Ch. 6 we can also estimate the error in the L''(0,T; F{h7j))- 
norms, r > 2, for non-dispersive schemes; for example for the classical three-point second 
order approximation of the laplace operator. In this case, in contrast with the good proper- 
ties of dispersive schemes, for iJ*(M)-initial data with small s, 1/2 — 1/r < s < 4-f- 1/2 — 1/r, 



(1.7) 
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the error losses a factor of order /i3/2(i/2-i/r) ^[^i^ respect to the case -L°°(0,T; PiJiL)) which 
can be handled by classical energy methods (see Example 1 in Section [3. 2p . Summarizing, we 
see that dispersive property of numerical schemes is needed to guarantee that the convergence 
rate of numerical solution is kept in the spaces L'^(0,T; /''(/iZ)). 

In the the context of the NSE we prove that the dispersive methods introduced in this 
paper converge to the solutions of NSE with the same order as in the linear problem. To be 
more precise, in Theorem 15.41 we prove a polynomial order of convergence, /i'^/^, in the case 
of a dispersive approximation scheme of order two for the laplace operator for initial data 
H'^iW^) when < s < 4. In the case of the classical non-dispersive schemes this convergence 
rate can only be guaranteed for smooth enough initial data, /7'*(M), 1/2 < s < 4 (see Theorem 

EI]). 

In Section[6]we show that non-dispersive numerical schemes with rough data behaves badly. 
Indeed, when using non-dispersive numerical schemes, combined with a i?^(]R)-approximation 
of the initial data ip G //*(R)\ff^(R), one gets an order of convergence | \ogh\~^/'^^~^^ which 
is much weeker than the /I'^/^-one that dispersive schemes ensure. 

The paper is organized as follows. In Section [2] we first obtain a quite general result which 
allows us to estimate the difference of two families of operators that admit Strichartz estimates. 
We then particularize it to operators acting on discrete spaces P{h'Z), obtaining results which 
will be used in the following sections to get the order of convergence for approximations of 
the NSE. In Section [3] and Section H] we revisit the dispersive schemes for LSE introduced in 
[ini [T71 [181 [20] which are based, respectively, on the use of artificial numerical viscosity and 
a two-grid preconditioning technique of the initial data. 

Section [5] is devoted to analyze approximations of the NSE based on the dispersive schemes 
analyzed in previous sections. Section [6] contains classical material on conservative schemes 
that we include here in order to emphasize the advantages of the dispersive methods. Finally, 
Section [7] contains some technical results used along the paper. 

The analysis in this paper can be extended to fully discrete dispersive schemes introduced 
and analyzed in [H] and to the multidimensional case. However, several technical aspects 
need to be dealt with carefully. In particular, one has to take care of the well-posedness of 
the NSE (see f5]). Furthermore, suitable versions of the technical harmonic analysis results 
employed in the paper (see, for instance. Section [7]) would also be needed (see ^3j). This will 
be the object of future work. 

Our methods use Fourier analysis techniques in an essential manner. Adapting this theory 
to numerical approximation schemes in non-regular meshes is by now a completely open 
subject. 

2. Estimates on linear semigroups 

In this section we will obtain LfL^. estimates for the difference of two semigroups SAit) and 
Ssit) which admit Strichartz estimates. Once this result is obtained in an abstract setting 
we particularize it to the discrete spaces /^(/iZ). 

2.1. An abstract result. First we state a well-known result by Keel and Tao [23] . 

Proposition 2.1. (\23\, Theorem 1.2) Let H he a Hilbert space, {X,dx) be a measure space 
and U{t) : H — )• L'^{X) be a one parameter family of mappings with i G M, which obey the 
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\mt)f\\LHX)<C\\f\\H 



energy estimate 
(2.1) 

and the decay estimate 

(2.2) WUmisygho^^x) < C\t - s|-"|bllLi(x) 
for some a > 0. Then 

(2.3) \\U{t)f\\L.^^^LriX))<C\\f\\H, 



(2.4) 



(2.5) 



(C/(s))*F(s,.))ds 



H 



[ U{t-s)F{s)ds 
Jo 



< C\\F\ 



Li{R,L^(X)) 



for all {q,r) and {q,f), a-admissible pairs. 



The following theorem provides the key estimate in obtaining the order of convergence 
when the LSE is approximated by a dispersive scheme. 

Theorem 2.1. Let {X,dx) be a measure space, A : D{A) L'^{X), B : D{B) L'^{X) 
two linear m-dissipative operators with D{A) D{B) continuously and satisfying AB = 
BA. Assume that {SA{t))t>o CLnd {SB{t))t>o the semigroups generated by A and B satisfy 
assumptions ()2.ip and ()2.2p with H = L^(X). Then for any two a-admissible pairs {q,r), 
{q, r) the following hold: 

i) There exists a positive constant C{q) such that 

(2.6) \\SA{t)ip - SB{t)Lp\\L<,{I,L-{X)) < CiQ)^^^{\MLHX),\I\\\i^ - B)ip\\L2(x)} 

for all bounded intervals I and ip G D{A) n D{B). 

ii) There exists a positive constant C{q,q) such that 
ft ft 



(2.7) 



SA{t - s)f{s)ds 



[ SB{t-s)f{s)ds 

Jo 



Li{I,L'-{X)) 



< C(g, g) min 1 1 1 / 1 1 , I /| 1 1 ( A - S) / 1 1 I 

for all bounded intervals I and f £ L^'{I, U' {X)) such that {A - B)f G L'?'(/, U'{X)). 

Proof of Theorem \2.1{ Using that the operators 5"^ and Sb verify hypotheses (12. ip and (12.2 
of Proposition O with H = L^iX), by ([O]) we obtain 



(2.8) 

and, by (1231) . 
(2.9) 



\\SA{t)ip - SB{t)ip\\Li{i,L-ixy) < C{<i)M\lhx) 



f SA{t-s)f{s)ds- f SB{t-s)f{s)di 

Jo Jo 



Li{R,L'-{X)) 



<Ciq,q)\\f\\ 



Li' (M, L?'(X))- 



In view of ()2.8|) and ()2.9p it is then sufficient to prove the following estimates: 
(2.10) \\SA{t)^ - SB{t)ip\\L.ii,Lrix)) < Ciq)\I\\\iA - BMl2^x) 
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and 
(2.11) 

/ SAit-s)fis)ds- [ SB{t-s)f{s)ds <C{q,q)\I\\\{A-B)f\\^,,^j^^,,^^^y 
In the case of (I2.10p we write the difference 5^(0 — Sb{-) as follows 

(2.12) SA{t)ip - SB{t)ip = [ SB{t-s){A-B)SA{s)ipds. 

Jo 

In order to justify this identity let us recall that for any (p G D(A) ^ D{B) we have that 
u{t) = SAit)ip £ C([0,oo),L»(A))nCi([0,(X)),L2(X)) and v{t) = SB{t)^ £ C{[0,oo), D{B))n 
C^([0, oo), L^(X)) verify the systems = Au,u{0) = ip, and vt = Bv,v{0) = ip respectively. 
Thus w = u-v £ C{[0,oo),D{B)) n C^{[0,oo), L^{X)) satisfy the system wt = Bw + {A - 
B)u,w{0) = 0. Since {A - B)u E C([0, oo), L2(x)) we obtain that w satisfies (f2l^ . 

Going back to (I2.12p and using that A and B commute we get the following identity which 
is the key of our estimates: 

(2.13) SA{t)^-SB{t)v= [ SB{t-s)SA{s){A-B)^ds. 

Jo 

We apply Proposition 12. 1 1 to the semigroup S'_b(-) and function F{s) = SAis){A — B)ip in this 
identity and, by (j2.5p with f = 2 and g = oo, we get 



(2.14) \\SA{t)ip-SB{tMLHI,Lr(X))<C{q)\\SA{s){A-BML^^J^r.HX)) 

<C{q)\I\\\{A-BML^xy 
Thus, (|2.10p is proved. As a consequence (|2.8p and (|2.10p give us (|2.6p . 

We now prove the inhomogenous estimate (j2.1ip . Using again (j2.13p we have 

SA{t-s)f{s)-SB{t-s)f{s)= f ' SB{t-s-a)SA{<y){A-B)f{s)da. 

Jo 

We integrate this identity in the s variable. Applying Fubini's theorem on the triangle {(s, o") 
0<s<t,0<(T<t — s} and using that A and B commute, we get: 

(2.15) A/(t) := /* SA{t - s)f{s)ds - f SB{t - s)f{s)ds 

Jo Jo 

pt pt—s 

= SB{t- s-a)SA{(r){A- B)f{s)dads 

Jo Jo 

= SB{t- s- a)SA{(T)ds{A- B)f{s)da 

Jo Jo 

= / SAia) / SB{t- s- a){A- B)f{s)dsda 
Jo Jo 



I SA{t-(j) I Sb{(7 - s){A- B)f{s)dsda 
Jo Jo 

f SA{t-c7)K^{A-B)f{a)da 
Jo 
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where 



Aig{t) = [ Seit - T)g{T)dT. 
Jo 



Applying the inhomogeneous estimate p.Sp to the operator Sa{-) with {q',r') = (1,2) we 
obtain 

(2.16) \\Afh.^i,Lrix))<CXq)\\Ai{A-B)fh^i,LHX))<C{q)\I\\\Ai{A-B)fh^ 
Using again ()2.5p for the semigroup S'_b(-), F = {A — B)f and (g, r) = (oo,2) we get 

(2.17) ||Ai(^ - B)fU^^j^LHx)) < emu - 

Combining and ([^TT]) we deduce ([^TT]) . Estimates (gS]) and (^TH} finish the proof. □ 

Remark 2.1. We point out that, in the the proof of the following estimate 

\\SAit)ip - SBitML.ilJr^x)) < C^('Z)|/|||(^ - BMl2^x), 

in view of (12.130 and (12.140 . we do not need that the two operators S'yi(t) and Ssit) admit 
Strichartz estimates. Indeed, it is sufficient to assume that only one of the involved operators 
admits Strichartz estimates and the other one to be stable in L^(X). 

2.2. Spaces and Notations. In this section we introduce the spaces we wih use along the 
paper. The computational mesh is hZ = {jh : j € for some h > and the P(/iZ) spaces 
are defined as follows: 

F{hZ) = {ip :hZ^C: \\(p\\iP(hZ) < ^} 

where 

{hY^HjhWy' l<p<oo, 

\ip{hz) = \ jez 

sup \u{jh)\ p = oo. 

On the Hilbert space P{hZ) we will consider the following scalar product 



{u,v)h =Re(h^u{jh)v{jh) ^ 

When necessary, to simplify the presentation, we will write {ipj)j^z instead of {ip{jh))j^z- 
For a discrete function {ip(jh)}ji^z we denote by (p its discrete Fourier transform: 

(2.18) ^(e) = /i5^e-^-«V(iM- 

iez 

For s > and 1 < p < oo, Ty*'P(]R) denotes the Sobolev space 

w''P{m) = {ipe s'{m) (/ - A)^/V e 

with the norm 

= ll((l + ICP)'*^^^)^IIlp(R) 

and by H'{R) the Hilbert space W'''^{R). 

The homogenous spaces W^'^{M), s > and 1 < p < c«, are given by 
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endowed with the semi-norm 

llv'llvy^.pCR) = II(I^I'''?)^IIlp(ir)- 

If p = 2 we denote H'{R) = W'^'^iR). 

We will also use the Besov spaces both in the continuous and the discrete framework. It is 
convenient to consider a function r]o G Cc(M) such that 

1 if lel<l, 
if |e|>2, 

and to define the sequence {'i]j)j>i S 5(M) by 



in order to define the Littlewood-Paley decomposition. For any j > we set the cut-off 
projectors, Pj(p, as follows 

(2.19) Pj^ = {rj.ipy. 

We point out that these projectors can be defined both for functions of continuous and discrete 
variables by means of the classical and the semi-discrete Fourier transform. 

Classical results on Fourier multipliers, namely Marcinkiewicz's multiplier theorem, (see 
Theorem I7.ip show the following uniform estimate on the projectors Pj: For all p £ (l,oo) 
there exists c{p) such that 

(2.20) \\PM\LPiR) < c{p)\\ip\\LPiR), yip e LP(R). 

We introduce the Besov spaces Bp2(M) ior 1 < p < oo hy 2 = {u € S'{R) : ||ti||_Ba ^(IR) ^ 
00} with 

^ 1/2 

lklb;,2(R) = \\Pou\\lp{K) + {^'^^'^\\PjU\\lp(R)) ■ 

Their discrete counterpart Bp2{h'^) with 1 < p < 00 and s G M is given by 



Pp,2ih^) = {u ■■ \\u\\B;^^{hZ) < 00}, 

with 

°° 1/2 

(2-21) = \\Pou\\iP(hz) + [^'^'^^'WPjuW'iPihZ)) ' 

where Pju given as in ()2.19p are now defined by means of the discrete Fourier transform of 
the discrete function u : /iZ — )• C. 

We will also adapt well known results from harmonic analysis to the discrete framework. 
We recall now a result which goes back to Plancherel and Polya [27] (see also [33], Theorem 17. 
p. 96, and the comments on p. 182). 

Lemma 2.1. ([27], p. 157) For any p G (l,oo) there exist two positive constants A{p) and 
B{p) such that the following holds for all functions f whose Fourier transform is supported 
on [— vr, vr]; 

(2.22) A{p) < [ Ifix^dx < B{p) ^ \f{m)\P. 

meZ m£Z 
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This result permits to show, by scahng, that, for all /i > 0, 
(2.23) < ||/|Up(k) < i?(p)^/P||/||/P(,z) 

holds for all functions / with their Fourier transform supported in [— vr/Zi, vr/Zi]. 

For the sake of completeness we state now the discrete version of the well known uniform 
L^-estimate (|2.20p for the cut-off projectors Pj. 

Lemma 2.2. For any p G (l,oo) there exists a positive constant c{p) such that 

(2-24) \\PjLp\\iv(^hZ) < c{p)\\(p\\iP(hZ) 

holds for all ip G F^hl?), j > 0, uniformly in h > 0. 

Proof. For a given discrete function ip we consider its interpolator (p defined as follows: 

J-Tv/h 

Thus, by ([2:23]) we obtain 

\\Pjf\\lP{hZ) < ciPM^jfTWlPm = c{p)\\Pj'P\\LP(R) < c{p)\\'p\\lp{R) < c{p)\\(p\\lP(hI.)- 

□ 

We recall the following lemma which is a consequence of the Paley-Littlewood decomposi- 
tion in the x variable and Minkowski's inequality in the time variable. 

Lemma 2.3. ([28\, Ch. 5, p. 113, Lemma 5.2) Letr] e C^{R) and Pj be defined as in (f2l9]) . 
Then 

(2.25) U\\h{R,Lr(R)) ^^\\Pj^\\Li(R,LriR)) if 2 < r < oo and 2 < q < CO 
and 

(2.26) \\PMh{R,Lr{R)) ^ II^IIL(R,L'-(R)) if'^<r<2andl<q<2 

j>0 

hold for all ijj £ L9(M, L'^(M)). 

Applying the above result and Lemma |2. II to functions with their Fourier transform sup- 
ported in [— 7r//i, 7r//i], as above, we can obtain a similar result in a discrete framework. 

Lemma 2.4. Let t] G C^iR) and Pj defined as in (IZTOD . Then 

(2-27) U\\li(R,ir^hZ)) ^ Yl \\PMLi(R,lr^hZ)) if2<r<ooand2<q<oo 

j>0 

and 

(2.28) Y ll^iV'llL(R,r(fez)) ^ M\h{R,r{hZ)) if I < r <2 and I < q <2 

hold for all G L''(]R, V{h'L)), uniformly in h > 0. 
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2.3. Operators on /P(/iZ)-spaces. In the following we apply the results of the previous sec- 
tion to the particular case X = KL. We consider operators j4/iwith symbol a^, ■ [— vr//i, vr/Zi] — )• 
C such that 

J-Tv/h 

Also we will consider the operator |V|* acting on discrete spaces P{h7j) whose symbol is given 

by ler- 

The numerical schemes we shall consider, associated to regular meshes, will enter in this 
frame by means of the Fourier representation formula of solutions. 

Theorem 2.2. Let Ah,Bh ■ l'^{hZ) — )• /^(/iZ) be two operators whose symbols are Oh and bh, 
ibh being a real function, such that the semigroups they generate, (S'A^(i))t>o o.'i^d {SB^{t))t>0! 
satisfy assumptions (12. ip and (12. 2p with some constant C , independent of h. Finally, assume 
that for some functions {^{k,h)}k^p , with F a finite set, the following holds for all ^ £ 
[—■K/h, -n/h]: 

(2.29) i«h(e)-feh(e)i<5]Mfc,/i)ici'. 

feeF 

For any s > 0, denoting 

(2.30) e(s,/i) = J]M^,/i)™'^WM}^ 

the following hold for all {q,r), {q,r), a-admissible pairs: 

a) There exists a positive constant C{q) such that 

(2.31) \\SAh{t)'f - SBSt)'f\\Li{ij-{hi)) < C(9)e(s,/i)max{l,|/|}||93||B| 2(fez) 

holds for all ip G i?! 2(^^) uniformly in h > 0. 

b) There exists a positive constant C{s,q,q) such that 

(2.32) f SA,{t-a)f{a)da- f S B,{t - a) f {a)da 

Jo Jo Li{I,l^(h1)) 

< C(s, q, q)e{s, h) max{l, |/|}||/|| Li'{I,Bp ^(hZ)) 

holds for all f £ L^'{I, Bp^^ihZ)). 

Remark 2.2. The assumption that the semigroups {SAf^{t))t>o o-i^d (S^^ (t))t>o, satisfy ()2.ip 
and (|2.2p with some constant C, independent of h, means that both of them are l'^{h'L)- stable 
with constants that are independent of h and that the corresponding numerical schemes are 
dispersive. 

Taking into account that both operators, and Bh, commute in view that they are asso- 
ciated to their symbols, the hypotheses of Theorem \2.1\ are fulfilled. They also commute with 
|V| and Pj which are also defined by a Fourier symbol. 

Assumption (12.29P on the operators Ah and Bh implies 

\\iAh-BhMi2^hZ) <^aik,h)\\\V\''ip\\i2^hZ)- 

keF 

However, this assumption is not sufficient to obtain a similar estimate in V [liL) -norms, r 7^ 2. 
As we will see this will be an inconvenient in obtaining ()2.32p as a consequence of ()2.7p . 
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The requirement that ibh is a real function is needed to assure that the semigroup generated 
by Bh, Sb^, satisfies 

Ss^it -cr) = SBjt)SB^i-cr) = SBjt)SB^i(j)*, 

identity which will be used in the proof. 

In Section\^we will give examples of operators Ah and verifying these hypotheses. In all 
our estimates we will choose bh(^) = , which is the symbol of the continuous Schrddinger 
semigroup. 

Proof of Theorem \2.iA We divide the proof in two steps corresponding to the proof of (|2.3ip 
and ()2.32p respectively. 

Step I. Proof of (lOTI) . We apply inequality ([2:25|) to the difference SA^{t)'f - SBi,{t)ip: 

1 /2 

WSAd^)^ - SBi,it)^p\\Li{I;lr(hZ)) < [^WPjSAd^)^- Pj 33^^)^1111 {1,1- (hZ))) ■ 

Using that Pj commutes with Sa^{-) and Sb,^{-) we get: 

(2.33) \\SA^{t)(p - SB,St)'P\\Li{I;lr(hZ)) < [^\\{SA^{t) - SB,St))Pj^\\h(Ijr(^f,Z))) ' ■ 

In order to evaluate each term in the right hand side of (|2.33|) we apply estimate (|2.6|) to 
the difference SA,^i-) — ^B^i') when acting on each projection Pjip. Thus, using hypothesis 
I^I7B\i we obtain: 

(2.34) \\SA,{t)Pj^ - SB,{t)PMLiii,nhZ)) 

< C{q)m&x{\I\,l}min{\\Pj(p\\i2^h^), \\{Ah - Bh)Pjf\\i2^hZ)} 

< C{q)max{\I\,l}mm^\\Pjip\\i2(^hZ),'^Kk,h)\\\V\''Pjip\\i2^hZ)} 

k&F 

< C{q) max{|/|, 1} ^ min | ||P,-(^||;2(ft2), /i(/c, /i)2-'^||Pj99||;2(;jZ)} 

fceF 

<C{q) max{|/|,l}||P,-99||i2(K)^min{l,^(A;,/i)2-''^}. 

fceF 

Going back to estimate (|2.33p we get 
\\SAu{i)^-SBu{i)^\\Limr{hX)) 

< C{q) max{|/|, 1} ( ^ \\PM\hm E {l' /^'(^' ''^ ■ 

We claim that for any j > the following holds 

(2.35) ^min{l,^2(^^;j)22-''^} < ^ ^(A:, /i)""'^^2s/fc,2}22is 
fceF k&F 

for all s > 0. 
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Assuming for the moment that the claim ()2.35p is correct we deduce that 
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\\SA,Si)'f-SBh{'t)'f\\Li{I;l-{hZ)) 

k(^F j>0 
k&F j>0 

< C{q,F)max{\I\,l}e{s,h)\\ip\\Bi^(R). 
We now prove (|2.35p by showing that 

(2.36) min{l,/i2^*^} < ^"^MVfc,i}2i^ 

holds for all // > and j > 1. It is obvious when fi > 1. It remains to prove it in the case 
/X < 1. For any |^| > 1 we have the following inequalities 

min{l,/i|e|^} < min{l,^|^|*^}™'^^"/^'^> = niin{l, 

^ ^^min{s/fc,l}|^|fcmin{s/fc,l} ^ ^min{s/fc,l} 

Applying this inequality to ^ = 2-', j > 0, we get ()2.36p and thus (j2.35p . The proof of the 
first step is now complete. 

Step II. Proof of (j2.32p . Let us denote by Ah the following operator: 

Ahf{t)= f SAdt-^)fi^)d^- f SBdt-^)fi^)d<r. 
Jo Jo 

As in the case of the homogenous estimate (j2.3ip . we use a Paley-Littlewood decomposition 
of the function /. Inequality (12.27P and the fact that Ah commutes with each projection Pj 
give us 

(2.37) WAhfWl^ijr^hZ)) < ^i'l)Y.\\^j^^^'f'>\\Li{i,i-{hZ)) = cil)^\\MPjf)\\li{i,ir(^hZ))- 

j>0 j>0 

We claim that each term A{Pjf) in the right hand side of (I2.37P satisfies: 
(2.38) 

\\^h{Pjf)\\Li{I,l^(hZ)) 



<c((?,g)max{l,|/|}min|||P,/||^,-.(,_^,.(^2;))'E/^(^'^)lll^l'^i/ll 

I k&F 



Li'{I,l^'{hZ)) ( ■ 



12 L. I. IGNAT AND E. ZUAZUA 

In view of (j2.36p . the above claim implies 

(2.39) ||A,(P,/)|| 

< c(g,g)max{l,|/|}mini 2^ /^(^' ^)2^''ll^i/llL4'a,«^'{hZ)) 

= c(g,g)max{l,|/|}||Pj/||^,-/(j ^min{l,^(A:,/i)2J''} 

k&F 

< c(g,glmax{l,|/|}||P,-/||^,,(,^,,,(,^„ j;^(A;,/.)--W^-'i}2^-^ 

< c{q, q) max{ 1 , | /| }e(s , /i)2^^ 1 1 / 1 1 ^,-/ ^^^-j) . 
Estimates (jOTl) and (jTM]) give us 

1 /2 

(2.40) \\Khf\\L.ii,irihZ)) < ciq,q)max{l,\I\}e{s,h)(^Y.2^^'\\Pjf\\l,,^j^^,,^^^^^) . 

j>0 

Using that q' < 2, we can use the reverse Minkowski's inequality in L^' /"^{I) to get 



^'^'^^"'W^i (I^lf {hi)) ~ X] ll^^^''ll^i-^llp--'(M))llL9'/2(/) ^ II ^'^'^^''\\^jf\\'ir'(h'E)) 

j>0 j>0 j>0 



L'5'/2(7) 



< 



j>0 



2 2 



II./ lli^'(J,B| .,(;,z))• 



By ([TIO]) we get 

l|Ah/llL9(/.''"(/*z)) - ^) max{l, |/|}e(s, /i)||/|| 

which finishes the proof. 

In the following we prove (j2.38p . Using that both operators Sa^ and Sb^ fulfill uniform 
Strichartz estimates, it is sufficient to prove that, under hypothesis (j2.29p . the following 
estimate holds for all functions / G L"? (/, /*" (/iZ)): 

(2.41) W^hfWmijr^hZ)) <c(g,g~)|I|^a(/e,/i)|||V|VllL^'(/,,*='(,,z))- 

fce-F 

We point out that, in general, this estimate is not a direct consequence of (j2.7p since, under 
assumption (j2.29p . we cannot establish the following inequality (of course, in the particular 
case f' = 2 this can be obtained by Plancherel's identity) 

ll(^, -p,)/lL,-' (,z))<E"(^'^)lll^l'/ll 

Identity (j2.15p gives us that 

Ah/(t) = r SA,{t - s)Aih{Ah - Bh)f{s)ds 

where 



Li (IJf (hZ))- 
k£F 







Jo 
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The inhomogeneous estimate (j2.5p with {q',r') = (1,2) shows that 
(2.42) W^hfhHu^ihZ)) < c{q) \\Aih{Ah - Bh)f\\mi,i2(^hZ)) ■ 

Using that satisfies 5*6^ (t — a) = S'B^(t)S'B^(— o") = SB^{t)SB^{o-)* and that it commutes 
with Ah we get 

Aih{Ah - Bh)f{t) = SB,St){Ah - Bh) f SBMTf{^)d<^- 

Jo 

Thus, using the uniform stability property, with respect to h, of the operators Sb^' 



hV)\\l'HhZ)~^l2{hZ) 



< 1 



and hypothesis (|2.29p we get 

(2.43) \\Aih(.Ah-Bh)f\\L^^i,mhZ)) < \\{Ah-Bh) [' SbA^T f(.^)d^ 

II Jo 



<y2a{k,h) |vi^' fsBMrf{^)d<T 

Using that B^ and |V| commute, estimate (j2.4p with [/(•) = Sb,^{-) gives us that 

\\Aih{Ah-Bh)f\\Lm,mhZ)) < \I\Y,a{k,h) f SBAsr\V\''f(.cT)da 

k&F -^0 

< |/| Va(fc,/i)sup / SBMTM''f{^)dT 

JCI J J 



L°°{I,P{hZ)) 



<c(g)|/|^a 

keF 

Thus, by (I2.42p we obtain (I2.4ip which finishes the proof. 



/2{/lZ) 



viVll 



L?' {I,r'{hZ)). 



□ 



3. Dispersive schemes for the linear Schrodinger equation 

In this section we obtain error estimates for the numerical approximations of the linear 
Schrodinger equation. We do this not only in the /^(/iZ)-norm but also in the auxiliary spaces 
that are needed in the analysis of the nonlinear Schrodinger equation. 

3.1. A general result. The numerical schemes we shall consider can all be written in the 
abstract form 

r iu'l{t) + AhU^ = 0, t>0, 

(3.1) <^ 

We assume that the operator Ah is an approximation of the 1 — d Laplacian. On the other 
hand, T/^c^ is an approximation of the initial data 93, Th being a map from L^(M) into P{hZ) 
defined as follows: 

(3.2) {Th^m) = e^^'^^md^- 

J-TT/h 
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Observe that this operator acts by truncating the continuous Fourier transform of on the 
interval (— vr/Zi, vr/Zi) and then considering the discrete inverse Fourier transform on the grid 
points KL. 

To estimate the error committed in the approximation of the LSE we assume that the 
operator A^, approximating the continuous Laplacian, has a symbol ah which satisfies 

(3.3) Mi) -e\<Y. «(^' ^ e [- J' I ' 

for a finite set of indexes F. As we shall see, different approximation schemes enter in this 
class for different sets F and orders k. 

This condition on the operator suffices to analyze the rate of convergence in the 
L°°{—T,T; I'^ih'L)) norm. However, one of our main objectives in this paper is to analyze this 
error in the auxiliary norms L'?(— T, T; F{h'Z)) which is necessary for addressing the NSE with 
rough initial data. More precisely, we need to identify classes of approximating operators A^ 
of the 1 — d Laplacian so that the semi-discrete semigroup exp{itA}i) maps uniformly, with 
respect to parameter h, P{hZ) into those spaces. 

In the following we consider operators Ah generating dispersive schemes which are PiJiL)- 
stable 

(3.4) II exp(itA/j)(^||p(^z) < C'llv^llp(hZ), Vt > 

and satisfy the uniform l^{hZ) — l°°{h'L) dispersive property: 

(J 

(3.5) II ey.Y>{itAh)'p\\i^{hZ) < -^^yh{hz), Vt > 0, 

for all /i > and for all (p G l^{hZ), where the above constant C is independent of h. We 
point out that (13.41) is the standard stability property while the second one, (13. Sp . holds only 
for well chosen numerical schemes. 

Applying Theorem 12.21 to the operator Bh whose symbol is —i^'^ and to iAh, Ah being the 
approximation of the Laplace operator with the symbol a/i(^), we obtain the following result. 

Theorem 3.1. Let s >0, Ah satisfying ()3.3p . ()3.4p . (j3.5p . and {q,r) and {q,f) be two 1/2- 
admissible pairs. Denoting 

(3.6) e{s, h) = J2 «(^' /i)"^^^^^^'!}, 

the following hold: 

a) There exists a positive constant C{q) such that 

(3.7) II exp(itA)T/,v3 - T/^ exp(ii9^)(/7||i5(o,r; r(hz)) < max{l, r}C(g)e(s, /i)||(/?||h=(]r) 
holds for all e iJ''(M), T > and /i > 0. 

h) There exists a positive constant C{q,q) such that 



(3.8) 



/ exp{i{t-a)Ah)Thf{a)da- [ T/, exp(i(t - c7)92)/(a)a!a 
Jo Jo 

<Ciq,q)T. 

holds for allT>0, f e L9'(0, T; Bf, ^(R)) and h > 



L'J(0,T;K(hZ)) 



C{q,q) max{l,r}e(s,/i)||/||^^/(o ,j,.^.^ 
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Remark 3.1. In the particular case when {q,r) = (c«,2) and the set F of indices k entering 
in the definition (I3.6P of e{s, h) is reduced to a simple element, the statements in this Theorem 
are proved in [30j (Theorem 10.1.2, p. 201): 

(3.9) II exp(zt^,,)T,,99 - T,,exp(iiO^)(^||i,x=(o,T;«2(M)) < C{l)Te{s,h)\\ip\\Hs{M.)- 

Remark 3.2. Observe that for s > sq = max{fc : k G F} the function s — )• e{s,h) is 
independent of the s-variable: 

e(s, k) = e(so, k) = a{k, h). 

k&F 

This means that imposing more than regularity on the initial data does not improve 

the order of convergence in (I3.7p and ()3.8p . 

Remark 3.3. In the case < s < sq, with sq as above, the estimate H^'^'{]^) — )• L°°(0, T; P{hZ)) 
in (|3.7p and the one given by the stability of the scheme — ?• L°°(0, T; /^(/iZ)), allow to 
obtain, using an interpolation argument, a weaker estimate: 

\\exp{itAh)Thip - Thexp{itdl)cp\\Lc^(^o^T;P{hZ)) < C{T)e{so,hy/"'\\Lp\\Hs(R). 

If the set F has an unique element then this estimate is equivalent to (|3.7p . However, the 
improved estimates ()3.7p and ()3.8p cannot be proved without using Paley-Littlewood's decom- 
position, as in the proof of Theorem \2.2l 

3.2. Examples of operators A^. In this section we will analyze various operators A^, which 
approximate the 1 — d Laplace operator d'^. 

Example 1. The 3-point conservative approximation. The simplest example of 
approximation scheme for the Laplace operator 9^ is given by the classical finite difference 
approximation A/j 

[6.10) [/^hU)j - . 

It satisfies hypothesis p.3p with F = {4} and a(4, /i) = h^. Thus, we are dealing with an 
approximation scheme of order two. Indeed, we have: 

vr VT" 
h'h. ' 

However, this operator does not satisfy (j3.5p with a constant C independent of the mesh size 
h, (see [IT], Theorem 1.1) and Theorem 13.11 cannot be applied. This means that we cannot 
obtain the same estimate as for second order dispersive schemes: 
(3.11) 

\\exp{itAh)Th(p - Thexp{itdl)(f\\L<i(o,T;r(hZ)) < C{q,T)\\(p\\Hs^^-^ 

However, using the ideas of Brenner on the order of convergence in the /^(/iZ)-norm, r > 2, 
([3], Ch. 6, Theorem 3.2, Theorem 3.3 and Ch.3, Corollary 5.1) we can get the following 



/)2 



sm 



</i'iei',vee 



r h^/^ SG (0,4), 
\ h^, s> 4. 



16 



L. I. IGNAT AND E. ZUAZUA 



estimates: 

52 



exp{itAh)Th(p-Th exp{itd^)(p\\L',(o,T; ir{hz)) 



<C{q,T)M 



<Ciq,T)\\^\\ ,,1 1 



h^, s>4 + l-f, 



/i^ s>4 + l-§ 



where we have used that i7^o(M) = ^aM^) -^r,oo(I^) when sq - 1/2 = s - l/r. 
Observe that in the case s G (0, 4) the above estimate guarantees that 

(3.12) II exp(itA/j)Tfe99-T/j exp(ita^)(/?||i,(o,T; i-{h-L)) 

<c(g,r)||^||^,,._i^^^/.i(^+^^)/.-l(^^). 

Moreover for any o" G (1/2 — 1/r, 4 + 1/2 — 1/r) we can find s G (0, 4) with cj = s + l/2— 1/r 
and using (j3.12p we obtain 

(3.13) \\exp{itAh)Thip - Thexp{itdlML,io,T-ir^hZ)) < C{q,T)\\ip\\H.^^)h^ . 

In the case of an approximation of order two one could expect the error in the above estimate 
to be of order /i*^/^ as in the L°^(0, T; P{h'Z)) case. But, here we get an extra factor of 
order /i~3/2(i/2-i/r) -^j^jch diverges unless r = 2, which corresponds to the classical energy 
estimate in L°°(0, T; L^(M)). This does not happen in the case of a second order dispersive 
approximation of the Schrodinger operator, where Theorem 13.11 give us an order of error as 

in dan]). 

Note that, according to Theorem 13. H this loss in the rate of convergence is due to the lack 
of dispersive properties of the scheme. 

Also we point out that to obtain an error of order h"^ in (f3T2]) we need to consider initial 
data in H^~^^~'^^'^ (M). So we need to impose an extra regularity condition of 1 — 2/r derivatives 
on the initial data if to assure the same order of convergence as the one in (13. lip for dispersive 
schemes. 

Example 2. Fourier filtering of the 3-point conservative approximation. Another 
example is given by the spectral filtering A/^^^ defined by: 

(3.14) Aft^^y? = A/i(l(_2ZL,2ZL)^)^, 7 < ^• 

In other words, Ah^-y is a discrete operator whose action is as follows: 

"^/^ 4 . 2^^h 



i.e. it has the symbol 



0'h,j{l;) - J^Sm [ — )i-(-yn/h,-y7T/h)- 
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In this case 



ahAO-n<c{l){ <c{j)h'e for all 



h' h 



Thus ,y constitutes an approximation of the Laplace operator A of order two and the 
semigroup generated by iA^^^ has uniform dispersive properties (see [13). Theorem 13.11 
which exploits the dispersive character of the numerical scheme, gives us 

h'/^, sG(0,4), 
h^, s > 4. 



\eyLp{itAh)Thip - TheyLp{itA)ip\\Lq^o^T;ir(hZ)) < C{q,T)\\ip\\Hs 



We note that using the same arguments based on /''(/iZ)-error estimates (given in [3j), as in 
the Example 1, we can obtain the same result only if r = 2 or assuming more regularity of 
the initial data if. 

This scheme, however, has a serious drawback to be implemented in nonlinear problems 
since it requires the Fourier filtering to be applied on the initial data and also on the nonlin- 
earity, which is computationally expensive. 

Example 3. Viscous approximation. To overcome the lack of uniform L"^ (I , F (hl^)) 

estimates, in |18) and [15] numerical schemes based in adding extra numerical viscosity have 
been introduced. The first possibility is to take = + ia{h)Ah with a{h) = /i^-i/^C^) 
and a{h) -^1/2 such that a{h) — )• 0. In this case p.3p is satisfied as follows: 

(3.15) 1^ sin^ (|^) + ^a{h)^ sin^ (|^) - e\ < + a(h)e. 

This numerical approximation of the Schrodinger semigroup has been used in [18] and [20] 
to construct convergent numerical schemes for the NSE. However, the special choice of the 
function a{h) that is required, shows that the error in the right hand side of (j3.15p goes to 
zero slower that any polynomial function of h and thus, at least theoretically, the convergence 
towards LSE, and, consequently to the NSE, will be very slow. Thus, we will not further 
analyze this scheme. 

Example 4. A higher order viscous approximation. A possibility to overcome the 
drawbacks of the previous scheme, associated to the different behavior of the /^(/iZ) — /°°(/iZ) 
decay rate of the solutions, is to choose higher order dissipative schemes as introduced in |15j : 

(3.16) Af, = Ah- ih^'^'^-^X-Ahr, m > 2. 
In this case, hypothesis p.3p reads: 

(3.17) 1^ sin2 {^) + 1) sin2 [^)y -e\< + /,2(--l)^2m_ 
Theorem 13.11 then guarantees that for any < s < 4 the following estimate holds: 

||exp(itAft)T;,(^-T;,exp(itA)y;||i,(o,T;r(hZ)) < max{l, r}(/i^/2 ^ /i("^-i)^/"^)||(^||H.(K) 

< max{l,r}/i"/2||(^||^,^j^^. 

Thus we obtain the same order of error as for the discrete Laplacian A^ = A^ but this time 
not only in the L°°{I; P(/iZ))-norm but in all the auxiliary L'^{I, r(/iZ))-norms. We thus get 
the same optimal results as for the other dispersive scheme in Example 2 based on Fourier 
filtering. 



18 



L. I. IGNAT AND E. ZUAZUA 



4. A TWO-GRID ALGORITHM 



In this section we analyze one further strategy introduced in [16], [18] to recover the unifor- 
mity of the dispersive properties. It is based on the two-grid algorithm that we now describe. 
We consider the standard conservative 3-point approximation of the laplacian: Ah = Ah. 
But, this time, in order to avoid the lack of dispersive properties associated with the high 
frequency components, the scheme will be restricted to the class of slowly oscillatory data 
obtained by a two-grid algorithm. The main advantage of this filtering method with respect 
to the Fourier one is that the filtering can be realized in the physical space. 

The method, inspired by |12j . is roughly as follows. We consider two meshes: the coarse 
one of size 4/i, h > 0, 4/iZ, and the finer one, the computational one, /iZ, of size h > 0. The 
method relies basically on solving the finite-difference semi-discretization on the fine mesh 
hTj, but only for slowly oscillating data, interpolated from the coarse grid 4/iZ. The 1/4 ratio 
between the two meshes is important to guarantee the dispersive properties of the method. 
This particular structure of the data cancels the pathology of the discrete symbol at the points 
±7r/2h. 

To be more precise we introduce the extension operator 11^'' which associates to any func- 
tion ip : 4hZ — 7- C a new function Tlf^ip : KL — t- C obtained by an interpolation process: 



(nf = (Pi,V)(j/i), i e 

is the piecewise linear interpolator of ifj. 
The semi-discrete method we propose is the following: 



where P^^ 



(4.1) 



iv!i{t) + Ahv!^ = i>0, 



The Fourier transform of the two-grid initial datum can be characterized as follows (see 
Lemma 5.2, [18]): 



(4.2) 



vr vr 
h'h 



where T/^hV^iO is the extension by periodicity of the function T^hf: initially defined on 
[— 7r/4/i, 7r/4/i], to the interval [— vr/Zi, vr/Zi], and 



(4.3) 



m(0 



e^^^ - 1 
4(e*« - 1) 



,p>2. 



The following result, proved in [16], guarantees that system (j4.ip is dispersive in the sense 
that the discrete version of the Strichartz inequalities hold, uniformly on /i > 0. 

Theorem 4.1. Let {q,r), {q,r) be two 1/2- admissible pairs. The following properties hold 
i) There exists a positive constant C{q) such that 

(4.4) l|e**^''ntVllLnM,r(/^z)) < c {q)\\n.f^ ^hz) 

uniformly on h > 0. 

a) There exists a positive constant C{d,r,r) such that 



(4.5) 



s<t 



<C(g,g)||nf/|| 



Li' 



,r'{hZ)) 
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for all f e L'?'(M, f{AhZ)), uniformly in h > 0. 

In the following lemma we estimate the error introduced by the two-grid algorithm. 

Theorem 4.2. Let s > and {q,r), {q,r) be two admissible pairs. 

a) There exists a positive constant C{q,s) such that 

(4.6) II exp{itAh)nl^T4hf-Thew{'i'tdl)(p\\Li{i;ir(hZ)) 

< C7(g,s)max{l,|/|}(/J-''^W2,2} ^;^min{.,l})||^||^^^^^^ 

holds for all f G H%R) and h > 0. 

b) There exists a positive constant C{q,q,s) such that 

(4.7) / eMi{t-s)Ah)Uj!'T4hf{s)ds- f Thexp{i{t - s)dl)f{s)ds 

Js<t Js<t 

< C(g,Q,s)max{l,|/|}(/i--W2,2} +;^min{.,l})||J||^,,^^^^^^ 



Remark 4.1. There are two error terms in the above estimates: /i™°i''/2'2} ^^^^ ^min{s,i}_ 
The first one comes from a second order numerical scheme generated by the approximation of 
the laplacian 5^ with Ah and the second one from the use of a two-grid interpolator. Observe 
that for initial data £ H^{M), s G (0,2) the results are the same as in the case of the 
second order schemes. Also, imposing more than H^{W) regularity on the initial data does 
not improve the order of convergence. This is a consequence of the fact that the two-grid 
interpolator appears. The multiplier m(^) defined in (|4.3p satisfies m{^) — as^~0 and 

then the following estimate, which occurs in the proof of Theorem 

\mm-i\''m)?di<{h\Mm^^^)f, 

-■K/Ah 

cannot be improved by imposing more regularity on the function (p. 



Proof of Lemma \4.2\ Case I. Proof of the homogenous estimate (j4.6p . Let us consider 
a'* acting on discrete functions as follows: 

J-w/h 

Note that A'* differs from the finite-difference approximation A^ on the fact that, in A'', ^"^ 
replaces the symbol 4//i^ sin^(,^/i/2) of A^. 

In view of the definition of A'^, we have 

exp{itA^)Th(p = Thexp{itdl)ip. 
Using the last identity, we write 

eMitAh)Ul^Th^ - Thexp{itdl)ip = exp(itAOnf T;,(^ - exp(itA^)T;,(^ 



hit) + hit) 

hit) = exp(^^Aft)^;^'^T4/^<^ - exp(itA^)nrT4/»(/J 



where: 
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and 

hit) = exp(itA'*)nf T4,^v' - exp{itA'')Thip. 
In the following we estimate each of them. 

Applying Theorem 12.21 to operators and A'* we get 

\\h\\LH0,T;irm) < /i"^^'^^^/2'=^Uax{l,r}||nf T4,^||B|,,(hZ) < /i"^"^^/'''>max{l,r}||^||^.(K). 

In the case of I2 we claim that for any s > 

(4-8) \\h\\Li{0,T;p-{hZ)) < ^™"^'*'"^^II'^IU''(M)- 

To prove this claim, we remark that the operator exp( itA^) satisfies (HH) and (g^]). Thus 
Proposition 12.11 guarantees that exp(ztA'^) has uniform Strichartz estimates and 

(4.9) \\h\\Li{o,T;r{hZ)) < \\^h"^'ihf - ThV\\i2(hz)- 
It is then sufficient to prove that 

(4.10) linf T4/,v^ - T^<^||p(;,^) < h"^'^^'''^ MlHsiu) 

holds for any s > 0. Actually it suffices to prove it for < s < 1. Also the cases s £ (0, 1) 
follow by intepolation between the cases s = and s = 1. We will consider now these two 
cases. 

The case s = easily follows since 
and 

\\T^h^\\l2(hZ) ^ IIV5||l2(R). 

We now prove (j4.10p in the case s = 1: 

(4.11) wnl^T^hv - Tfc(^llp(hz) < hiMmm- 

Using that 

it is sufficient to prove the following estimate 

(4.12) ||n^''T4,j(/9 - T4/i(/?||p(/j2) < 
The representation formula (j4.2p gives us that 

/7r/4/i 
\Ta(hi)-\\^\(p(i)\^di 

+ / \m{hO\^\T^(0?di. 

J-K/ih<\i\<Tx/h 

Using that |m(^) — 1| < |C| for ^ G [— 7r/4, 7r/4] we obtain 

/TV /Ah 
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Previous results on the Fourier analysis of the two- grid method (see [19], Appendix B) and 
the periodicity with period tt /2h of the function T/^}^'^{^) give us that 

4 



\m{hi)\''\T^HV{0?di 

7r/4ft.<|5|<7r//i 







1 


/ 

/7r/4/i.<|g|<7r//i 


4(e««^ - 


1) 



|T4h¥^(0l 



2 



J7r/4/i<|5|<7r//i J-7r/4h<g<7r/4/i 
J -■KlAh<e<-K/Ah 



-7r/4fe<g<7r/4/i 

We obtain that ()4.12p holds and, consequently, (j4.1ip too. Thus (|4.8p is satisfied for any 
positive s. 

Observe that the main term in the right hand side of (|4.13p is given by (|4.14p . and this 
estimate cannot be improved by imposing more than H^{U) smoothness on ip. 

Case II. Proof of the inhomogeneous estimate (|4.7p . We proceed as in the previous 
case by splitting the difference we want to evaluate as 

/ exp(i(t - s)^h)^fTAhf{s)ds - [ Th exp(i(t - s)dl)f{s)ds = h + h 

J s<t J s<t 

where 

;exp(i(t - s)Ah) - exp(i(t - s)A^'))nfT^hf{s)ds, 



s<t 



and 



l2= [ eMi{t - s)A''){UfT^hf(.s) -Thf{s))ds. 

Js<t 

In the case of Ji, applying Theorem 12.21 to operators Ah and A'', we get 

\\h\\Li{o,T-r{hi.)) < /i™"'^*/2,2} niax{l,T}||n^''T4/j/||^,-/(Q^^.g,^ ^(feZ))- 
Applying Theorem 17. II below to the multiplier m given by (j4.3p . for any s > we obtain that 

\\^h''^ihf\\Li'{0,T;Bp^{hZ)) - II Hi^' (0,T; B^, ^W) 

and then Ii satisfies: 

(4.15) ||/i||L.(o,T;r'-(/.z)) < /i°'''^^^/2'2Uax{l,r}||/||^,-,(o^^^^.^^(^)). 

In the case of I2 we claim that 

(4.16) l|/2||Lno,T;r(hZ))</i°"°^^''^l 



'2||L9(0,T;r{/iZ)) ^ \\J IIl9'(0,T;B-, ^(M))- 

To prove this claim we consider the cases s = and s = 1. When s £ (0, 1) we use interpolation 
between the previous ones. Also the case s > 1 follows by using the embedding Bp 2(1^) ^ 

The case s = follows from Proposition 12.11 applied to the operators Uh{t) = T/j exp(it9^). 
We now consider the case s = 1. Using Strichartz estimates given by Proposition 12.11 to 
the operator exp(itA'^) we get: 

l|-^2||L9{0,r;r(/iZ)) < W^h'^^Ahf - T^hfWLij' ^o^T;lf' (hZ))- 
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Theorem 17.11 applied to the multipher m gives us 

||n^'^T4/i/ - T^4hf\\Lq'(^o^T;lf'{hZ)) — / II L?' {0,T; fii, ^(R)) 

and 

IIT4/1/ - T^hf\\Lq'(^Q^T;lf'{hZ)) - ^II/IIl'?'(0,T;B1, ^(R))- 

Thus (|4.16|) holds for s = 1, and in view of the above comments, for all s > 0. 
Putting together (j4.15p and (|4.15p we obtain the inhomogeneous estimate (|4.7p . 
The proof is now complete. □ 

5. Convergence of the dispersive method for the NSE 

In this section we introduce numerical schemes for the NSE based on dispersive approxi- 
mations of the LSE. We first present some classical results on well-posedness and regularity 
of solutions of the NSE. Secondly we obtain the order of convergence for the approximations 
of the NSE described above. 

5.1. Classical facts on NSE. We consider the NSE with nonlinearity f{u) = \u\Pu and 
93 G H^{M). We are interested in the case of H^(R) initial data with s < 1. The following 
well-posedness result is known. 

Theorem 5.1. Let f{u) = \u\Pu with p G (0,4). Then 

i) (Global existence and uniqueness, [5j, Th. 4-6.1, Ch. 4, P- 109) 
For any Lp G L^(M), there exists a unique global solution u of (|1.2p in the class 

uG C(R,L2(R))nLf„^(M,L^(M)) 

for all 1/2- admissible pairs {q,r) such that 

\\u{t)\\L2(R) = \\(p\\l2{r), G M. 

a) (Stability, [5], Th. 4-6.1, Ch. 4, p. 109) Let ip and ip be two functions, and u and 

V the corresponding solutions of the NSE. Then for any T > there exists a positive constant 
C{T, \\'P>\\L2(R),\\tp\\L^(R)) such that the following holds 

(5-1) ll^i - v\\l^{0,T;L^R)) < C{T, MIl^r), IIV'IIl2(r))II95 - V'IIl2{R) 

Hi) (Regularity) Moreover if if £ H^{W), s G (0,1/2) then Theorem 5.1.1, Ch. 5, 
p. 147) 

u G C(M,-H''(M)) n Lf„^(M,5^"_2(^)) 
for every admissible pairs {q, r) . 

Also ifipe H'^{R) then u G C{R,H'^(R)) ([5\, Theorem 5.2.1, Ch. 5, p. I49). 

Remark 5.1. The embedding B'^^^i^) ^ W'^''{M.),r > 2, (see [5\, Remark I.4.3, p. I4) 

guarantees that, in particular, u G L^^^(M, Ty*''~(M)). Moreover, f{u) G L^^J^, B^, ^i^)) and 
for any < s < 1 (see [5], formula (4.9.20), p. 128) 

4-p(l-2s) 

(5-2) II/('")IIl9'{7,B», ^(K)) ^ l-^l II^IIl9(/,B^_2W)' 

The fixed point argument used to prove the existence and uniqueness result in Theorem 
15.11 gives us also quantitative information of the solutions of NSE in terms of the L^(M)-norm 
of the initial data. The following holds: 
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Lemma 5.1. Let (p G L^(M) and u he the solution of the NSE with initial data ip and 
nonlinearity f{u) = \uY'u, p G (0,4), as in Theorem \ 5.1[ There exists c{p) > and Tq = 
c(p)||(^||^2(]K)^ such that for any 1/2- admissible pairs {q,r), there exists a positive constant 
C{p, q) such that 

(5.3) \\u\\li{I;L-(R)) <C{p,q)\\ip\\L2(R) 

holds for all intervals I with \I\ <Tq. 

Proof of Lemma \5.1\ Let us fix an admissible pair {q,r). The fixed point argument used in 
the proof of Theorem 15.11 (see ([4], Th. 5.5.1, p. 15) gives us the existence of a time Tq, 

_ 4p 

To = C{p)\\(p\\^2[^y 

such that 

ll'"llL<'(o,ro;L'-(K)) < C{p,q)\\ip\\L2(^y 

The same argument apphed to the interval [{k — l)To,A;To], k > 1, and the conservation of 
the L^(M)-norm of the solution u of the NSE gives us that 

\H\Liak~l)To,kTo;L-{M.)) < C{p,q)\\u{{k - 1)Tq)\\l2(^k) = C(P, II</'IIl2(k) . 
This proves (j5.3p and finishes the proof of Lemma 15.11 □ 

5.2. Approximation of the NSE by dispersive numerical schemes. In this section 
we consider a numerical scheme for the NSE based on approximations of the LSE that has 
uniform dispersive properties of Strichartz type. Examples of such schemes have been given 
in Section [3] and Section HI 

To be more precise, we deal with the following numerical schemes: 

• Consider 

r iu1 + AhU^ = fiu''), t>0, 



(5.4) 



where A/^ is an approximation of A such that exp(itA/j) has uniform dispersive prop- 
erties of Strichartz type. We also assume that Ah satisfies Re(Mft,</?, ip)h < 0, K being 
the real part, and has a symbol a/i(^) which verifies 



(5.5) \ahiO-f\<^a{k,hm^U 



vr TT 
"h'h 



(5.6) 



• The two-grid scheme. The two-grid scheme can be adapted to the nonlinear frame as 
follows. Consider the equation 

in','' + A.nO''^ = nf /((nf )*nO''^), t > 0, 

n°''^(o) = nf 

where (nf )* : P{hZ) f{Ah'L) is the adjoint of : P{AhZ) P{hZ) and if'' is 
an approximation of f. 
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By [16], Theorem 4.1, for any p £ (0,4) there exists of a positive time Tq = 
TqUMlHr)) and a unique solution u^'^ G C(0,ro; P{hZ'^)) D Li{0,To; P+^ihZ'^)), 
q = 4(p + 2)/p, of the system (15. 6p . Moreover, u^'^ satisfies 

and 

(5-8) lk''llL9(o,To;«P+2(hZrf)) < c{TQ)\\li^ip^\\i2(^t^idY 

where the above constant is independent of h. 

With To obtained above, for any A; > 1 we consider u^'^ : [/cTq, {k + 1)Tq] — >• C the 
solution of the following system 



(5.9) 



.^'^ + A.n^''^ = nf /((nf )*n^''^), t G [kn, (k + i)ro]. 



u^^^ikTo) = Ufu^-^^^ikT^). 

Once, vf''^ are computed the approximation of NSE is defined as 

(5.10) u^{t) = n^''^(t), t G [kTo, {k + l)ro). 

We point out that systems (|5.6p and (j5.9p have always a global solution in the class 
C(M, P{hZ)) (use the embedding P{hZ) C l°°{hZ), a classical fix point argument and 
the conservation of the /^(/iZ)-norm). However, estimates in the L'^{0,T; r{h'L))- 
norm, uniformly with respect to the mesh-size parameter /i > 0, cannot be proved 
without using Strichartz estimates given by Theorem 14.11 Thus we need to take 
initial data obtained through a two-grid process. Since the two-grid class of functions 
is not invariant under the flow of system (15. 6p we need to update the solution at some 
time-step Tq which depends only on L'^(M)-norm of the initial data Lp. 

The following theorems give us the existence and uniqueness of solutions for the above 
systems as well as quantitative dispersive estimates of solutions v!^, similar to those obtained 
in Lemma |5. II for the continuous NSE, uniformly on the mesh-size parameter /i > 0. 



Theorem 5.2. Let p G (0,4), f{u) = \u\Pu and he such that Re{iAfi(p,f)h < and ()3.5p 
holds. Then for every (p^ G l'^{hZ), there exists a unique global solution G C(M, l'^{h'L)) of 
\5.4^ which satisfies 

(5-11) \\u^\\L^{R,l'2(h1)) < \W^\\l'^(h1)- 

Moreover, there exist c{p) > and C{p,q) > such that for any finite interval I with 

\i\<n = cip)y%lli;-^^ 

(5-12) \\u^\\Liii,r{hZ)) < C{p,q)\\'P^\\i2(^hZ), 

where {q,r) is a 1/2- admissible pair and the above constant is independent of h. 

Proof. Condition Re(i^/i</j, (/?)/i < implies the l'^{h'L) stability property (13. 4p . Then local 
existence is obtained by using Strichartz estimates given by Proposition 12.11 applied to the 
operator exp(it^/j) and a classical fix point argument in a suitable Banach space (see p8l and 
[20] for more details). The global existence of solutions and estimate (|5.1ip are guaranteed 
by the property ^e{iAhtp,ip)h < 0, and that Ke{if{uh),Uh)h = and the energy identity: 

(5.13) j^h\t)\\pihZ) = 2Re{iAhu\u% + 2Re{if{u''),u% < 0. 
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Once the global existence is proved, estimate ()5.12p is obtained in a similar manner as Lemma 
15.11 and we will omit its proof. □ 

Theorem 5.3. Let p G (0,4) and q = 4(p + 2) /p. Then for all h > and for every 
ip^ G P{4hZ), there exists a unique global solution G C(IR, l'^{hZ)) n lP+'^{hZ'^)) of 

()5.6p - ()5.10p which satisfies 

(5-14) ll'"''llL°°(M,/2(fcZ)) < W'^T'^^Wl'^ihl)- 

Moreover, there exist c{p) > and C{p,q) > such that for any finite interval I with 

\i\<n = cip)y%ll^^-^^ 

(5.15) h'^Wmup+^hZ)) < c(p,9)l|ntV|b2(;,z), 

where {q,r) is a 1/2- admissible pair and the above constant is independent of h. 

Proof. The existence in the interval (0, Tq), Tq = To(||9?'^||p(/j2)) for system (j5.4p is obtained 
by using the Strichartz estimates given by Theorem 14.11 and a classical fix point argument in 
a suitable Banach space (see [18] and [20] for more details). 

For any k > 1 the same arguments guarantee the local existence for systems (j5.9p . To 
prove that each system has solutions on an interval of length Tq we have to prove a priori that 
the /^(/iZ)-norm of does not increase. The particular approximation we have introduced 
of the nonlinear term in (j5.6p -( f5r9]l gives us (after multiplying these equations by vf''^ and 
taking the /^(/iZ)-norm) that for any t G [kT^, [k + 1)Tq] 

iW^'HrniHtz) = \\u'''''{kTo)\\i2^hZ) < \W''-''HkTo)\\mhZ) 

and then 

h''\t)\\iHhz) < \\u''\0)\\i2^hz) = \\uM\i2^hzy 

This proves (15.14p and the fact that for any k > 1 system (j5.9[) has a solution on the whole 
interval [kTo, [k + 1)Tq\. Estimate (|5.15p is obtained locally on each interval [/cTq, {k + l)ro] 
together with the local existence result. □ 



Let us consider the solution of the semidiscrete problem ()5.4p and u of the continuous 
one ()1.2p . In the following theorem we evaluate the difference between and ThU. 

Theorem 5.4. Let p G (0,4), s G (0,1/2), f{u) = \u\'Pu and be as in Theorem\5^ satis- 
fying dSSI). For any ip G H''{M.), we consider and u G L°°(M, iJ*(M)) n L^^"^(M, ^^+2,2(1^)); 
go = 4(p + 2)/p solutions of problems (j5.4p and (|1.2p . respectively. Then for any T > there 
exists a positive constant C{T, ||v'||l2(]R)) such that 

(5.16) 

\\u^ - T/iM||i9o(0,T;/P+2(/iZ)) + - T/jn||ioo(o^T;«2(M)) 

< \\^p\\l^r),p) e(s, h)\\u\\Lr=o^o,T;H={R)) + {h" + e(s, h)^ lhll?X\o,T;S;^2 
holds for all h > 0. 



2\ 



In the case of the two-grid method, the solution of system (15. 6p approximates the 
solution u of the NSE ()1.2p and the error committed is given by the following theorem. 
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Theorem 5.5. Let p ^ (0)4), s G (0,1/2), /(n) = \u\Pu. For any ip £ H^{M.), we consider 



andue L°°(M, i?^(M))nLf°^(M, 5^ , 2 2 W); Qo = 4(p+2)/p, solutions of problems (|5:6D - (|5lO]l 



ond (|1.2p . respectively. Then for any T > i/iere exists a positive constant C{T, WfWi'^iR)) 
such that 

(5.17) ll'u'' - Thu\\Lqo{0,T;lP+2{hZ)) + " T^hu\\L°" {0,T;P{hZ)) 

< C(r, ||(/?||i2(K),p) /l'*/^||'u||ioo(o,T;//''{R)) + + ll^llLt)\o,r;S;+2 2W) 

/loZds /or aZZ h > 0. 

Remark 5.2. Using classical results on the solutions of the NSE (see for example [4|, The- 
orem 5.1.1, Ch. 5, p. 147) we can state the above result in a more compact way: For any 
T > there exists a positive constant C{T, ||¥'||h=(r)) such that 

(5.18) \\u'^ - T/in||i9o(o,T;/p+2(hZ)) + 11^'' - T/in||ioo(o^T;/2{hz)) < C{T, \\(p\\H''{R))h^^'^ 
holds for all h > 0. 

Theorem 15.41 shows that if h^ < e{s, h) then the error committed to approximate the 
nonlinear problem is the same as for the hnear problem with the same initial data. As we 
proved in Section [3.21 for the higher order dissipative scheme = A/j — i/i^("*~^)(— A/j)*", 
m > 2, and for the two-grid method, £{s,h) = h^/'^ > h^. So these schemes enter in this 
framework. It is also remarkable that the use of dispersive schemes allows to prove the 
convergence for the NSE and to obtain the convergence rate for H'^(M) initial data with 
< s < 1/2. We point out that the energy method does not provide any error estimate in 
this case, the minimal smoothing required for the energy method being H'^(M), with s > 1/2 
(see Section [6] for all the details). 

In the following we prove Theorem 15. 4[ the proof of Theorem 15.51 being similar since 
the estimates in any interval (0, T) are obtained reiterating the argument in each interval 
{kTo, {k + 1)Tq), k >0, for some Tq = To{\\if\\ ^2 f^-^) in view of the structure of the scheme. 



Proof of Theorem \5.4[ The idea of the proof is that there exists a time Ti depending on the 
L^(M)-norm of the initial data: 

Ti:.min{l,||^||-^J'^f-^)}, 

such that the error in the approximation of the nonlinear problem 

evvh{t) = u^{t) - Thu{t), 

when considered in the -L'"'(0, Ti; ZP+^(/iZ))nL°°(0, Ti; /^(/iZ))-norm is controlled by the error 
produced in the linear part 

err';,*"(t) = exp(iM;,)T?,(/. - T;, exp(ita2)(^. 

In the following we denote by (g, r) one of the admissible pairs (00, 2) or {qo^p + 2). We 
now write the two solutions in the semigroup formulation given by systems (|5.4p and (jl.2p : 

u''{t) = eMitAh)T^h^ + i [ exp(i(i - s)Ah)f{u''{s))ds, 

Jo 
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and, respectively, 

Thu{t) = Th exp{itdl)ip + i [ Tt exp(i(t - s)d^)f{u{s))ds. 

Jo 

Thus 

(5-19) l|err/j||i9(o,T;r('i^)) - \\^^^''h'\\Li{o,T;r{hZ)) + l|err^''"'||L'j(o,T;r(hZ)) 

where, by definition, 

errr"(t) = T exp(i(t - s) A)/(n'^(s))ds - [\hexp{i{t - s)dl)f{uis))ds. 
Jo Jo 

For the linear part the error is estimated in Theorem 13.11 
(5.20) l|err|,™||i,(o,T;r{/iZ)) < C{q)e{s,h)mayi{T,l}\\ip\\H^^R). 

In the following we will estimate err^°". We write err;j°"(t) = I^{t) + I^{t) where 

I','{t) = I" exp(i(t - s)Ah) [f{u\s)) - ThfHs))) ds 



and 



Isit) = (exp(^(t - s)Ah)Thf{u{s)) - T;,exp(i(t - s)dl)f{u{s)))ds. 



Step I. Estimate of /g . For the last term, the inhomogeneous estimate (I3.8p in Theorem 
13.11 and estimate ()5.2p give us that 

(5.21) ||/3^*)llL^{o,T;r(^2;)) <C(g)e(5,/i)max{l,r}||/(tx)||^,,^^^^.^^ ^^^^ 



{p+2)',2 

< C{q)e{s, h) max{ 1, TjT^^^T^ II«IIlI(o,T;B: 



Step II. Estimate of We now prove the existence of a time Tq such that for all 
T <To, satisfies 

(5.22) \\l2{t)\\Li{0,T;l^{hZ)) 

< C{p)T'^-i\\eTih\\Lm{o,T-jp+HhZ))\Ml2(^R) + h'T^~^\\u\\^^^ 



P+2,2V 



The inhomogeneous Strichartz's estimate ()2.5p applied to the operators {exp{itAh))t>o 
shows that 
(5.23) 

\\llmLH0,T;nhZ)) < -T,/(n)||^,,^^^^^^(^^,,,^^^^^ 

< Cmf in') - /(T.n)ll^,, + CmfiT.u) - T./(.)||,,i(,^^^,,,,, 

We evaluate each term in the right hand side of ()2.3ip . In the case of the first one, applying 
Holder's inequality in time we get 

\\fiu')-f{T,n) 



L''0(0,T;Z(P+2)'(/iZ)) 

L90{0,T;«P+2(/iZ)) [\\'^"\\'LiO{0,T;lP+2{hZ)) + H "'-''J^IIl'JO (0,T; «P+2(/iZ)) ) ' 
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Let us now set Tq as it is given by Lemma |5. II and Theorem 

4p 

To ~ M\L2X^y 

Thus, by Theorem l5.1l Lemma [5.1l and Theorem l5.3l both and T/jU have their L'^{0, T; V'{h'L))- 
norm controlled by the L^-norm of the initial data: 

ll^''llLW(o,To;iP+2(/iZ)) < C(p)II</'IIl2(R) 

and 

I|T/i'U.||l<jo(0,To;Zp+2(/iZ)) < C'(p) II ^11 LTO (O.Tq; Lp+2(]R)) < C'(p) II 9^11 L2 (u) . 

These estimates show that for any T < Tq the following holds: 

(5.24) \\f{u^) - /(Th^^)ll^,^,(o,^^,(,+2/(,^)) < C(J,)T'~^u'' - T^n||^,o(o,T; ;p+2(,^)) ||v^||^,(j,). 

It remains to estimate the second term in the right hand side of ()5.23p . We will use now 
the following result which will be proved in Section [71 

Lemma 5.2. Let s G [0, 1], p > and f{u) = |n|^n. Then there exists a positive constant 
c{p, s) such that 

(5.25) ||/(T;,n)-T;,/(n)|| 

i(p+2y(^hx) — c(p, ■s)/i''||n||^^"^p_,.2j.jj-j 

holds for all u £ Ty^'P+2(R) andh>0. 

Using this lemma, Holder inequality in time and the embedding Bp^2 2i^) H^'''^^^(M) 
([5], Remark 1.4.3) we obtain: 

(5.26) ||/(T,^) - T./(^)||^,.(„^^^^,,,,,(^^^^ < c(p,.)/.^ri-f ||nr4„^^^^.„,,(^)) 

< c{p,s)h^T^ * II^IIl'?0(0,T;B»^2,2W)' 

Both (fOD and ([OGjl show that hit) satisfies ([5:22]) . 

Step III. Estimate of err/,. Collecting estimates ([OO]) . ([OT]) and ([5^22]) for both 
((7,r) = ((70)P + 2) and ((7,r) = (oo,2) we obtain that for any T < Tq the error err/, satisfies: 

(5.27) ||err/j||iqo(o,T;«p+2(/tZ)) + l|err/j||ioo(o,T;«2(/iZ)) 

<C{p) max{l,r}e(s,/i)||(^||H-(K) + C(p)||err/j||iTO{o,r;«p+2(/iZ))?'^"^'''^ll¥'lli2(K) 

+ h''T^ * II"IIl«)(0,T;B«^2_2(R)) + ^) "^^^"L^'^}^ " II'"IIl«)(0,T;B^^2,2W)' 

Now, let us set Ti < min{l,To} such that 

c{p)Tl-^/'Ml,^^^ < \. 

Then the error term err^ in the right hand side of ([5.27p is absorbed in the left hand side: 

l|err/i||L9o(o,Ti;ZP+2(hZ)) + l|errft||L°o(o,Ti;«2{/iZ)) 

< C(j))e{s,h)\\ip\\iisi^^-) + C{:p)\\u\f^^^^^^^ 
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We now obtain the same estimate in any interval (0, T). Using that the L^(R)-norm of the 
solution u is conserved in time we can apply the same argument in the interval [kTi, (k+l)Ti]: 

l|err/i||j;^90{fcri,{fc+l)Ti;/P+2(hZ)) + \\^''^''^h\\L°°{kTi,{k+l)Ti;l'^{hZ)) 

Let us choose T > and > 1 an integer such that (A^ — l)Ti < T < NTi. Thus 

l|err/i||2,9o(o,T;/P+2(feZ)) + \\^^^h\\L'=<'{0,T;l2(hZ)) 
N-l 

- X/ {\\'^^^h\\LiO{kTi,{k+l)Ti;lP+^{hZ)) + \\^^^h\\L°°{kTi,{k+l)Ti;P{hZ))^ 



k=0 

N-l N-l 

p+1 

L'J0(fcTi,(fc+l)ri;iJ;+2 2( 

fc=0 k=0 



< C{p)e{s,h) ^ \\u{kTi)\\Hs^u) + C{p)[h^+e{s,h)) ^ \\u 



Using that {p+ l)/qQ < 1 we have by the discrete Holder's inequality that 

ll"llL'JO(fcTi,(fc+l)Ti;B^ 2 2(R)) - II "IIlto (0,T; 2 



k=0 

Thus the error satisfies: 



|err/i||igo(0,T;/P+2(/iZ)) + \\'i'^'^h\\L'=°{0,T;P(hZ)) 

( \ i_E±l 

< iVe(s,/i)||n||Loo(o,T;H=(R)) + ( /i'' « II 

^llL'?0(0,T; _B 

< N 



e(s, h)\\u\\Lc^^Q^T;H''m + + ^)) 11^11^+^ 



i«0(0,T;B;+2, 



2^ 



< C{T,\\ip\\L2(^)) e(s,/l)||'u||ioo(o,T;//-(K)) + (^'^ + ^(S' ^)) ll^llLto\o,T; B;_^2 2(K 

This finishes the proof of Theorem 15.41 □ 

6. NONDISPERSIVE METHODS 

In this section we will consider a numerical scheme for which the operator has no uniform 
(with respect to the mesh size h) dispersive properties of Strichartz type. Accordingly we may 
not use L'^L'^ estimates for the linear semigroup exp{itAh) and all the possible convergence 
estimates need to be based on the fact that the solution u of the continuous problem is 
uniformly bounded in space and time: u £ L°°{{0,T); L°° (M.)). Thus, the only estimates 
we can use are those that the L^-theory may yield. When working with //''(R)-data with 
s > 1/2, using L°° (M; H'^ (M)) estimates on solutions and Sobolev's embedding we can get 
L^-estimates. 

There is a classical argument that works whenever the nonlinearity / satisfies 

(6.1) \f(^u)-fiv)\<Ci\u\P + \v\P)\u-v\. 

Standard error estimates (see Theorem 13.11 with the particular case {q,r) = (oo,2) or [30] . 
Theorem 10.1.2, p. 201) and Gronwall's inequality yield when < t < T: 

(6.2) \\u\t)-Thuit)\\i2^i,^) 

< h^/'^C{T){\\ip\\Hi(R) + lhllL"t\o,r:Hl(R)))«^^P(^ll'"llL-(0,T;Hi(M)))' 
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for the conservative semi-discrete finite-difference scheme. For the sake of completness we will 
prove this estimate in Section 16.11 

We emphasize that in order to obtain estimate ()6.2D we need to use that the solution u, 
which we want to approximate, belongs to the space L°°(M), condition which is guaranteed 
by assuming that the initial data is smooth enough. However, obviously, in general, solutions 
of the NSE do not belong to L°°(M) and therefore these estimates can not be applied. One 
can overcome this drawback assuming that the initial data belong to ff^(R) or even to //*(M) 
with s > 1/2 since in this case /f*(]R) Using -energy estimates and Sobolev's 

embedding we can deduce L°°-bounds on solutions allowing to apply (]6.2p . We emphasize 
that this standard approach fails to provide any error estimate for initial data in H'^(M) with 
s < 1/2. 

However, this type of error estimate can also be used for //*(M)-initial data with s < 1/2 
(or even for L^(M)-initial data), by a density argument. Indeed, given ip G H^(R) with 
< s < 1/2, for any 6 > we may choose ifs S H^CK) such that 

- ^5\\h={R) < ^■ 

Let us be the solution of NSE corresponding to ips- Obviously, ips being i?i(M)-smooth, 
we can apply standard results as (16. 2|) to us. On the other hand, stability results for NSE 
allow us to prove the proximity of u and us in (R) . This allows showing the convergence of 
numerical approximations of us, that we may denote by us^h, towards the solution u associated 
to ip as both (5—7-0 and /i — t- 0. But for this to be true h needs to be exponentially small 
of the order of exp(— 1/5) which is much smaller than the typical mesh-size needed to apply 
the results of the previous sections on dispersive schemes that required h to be of the order 
of 5^/'. 



6.1. A classical argument for smooth initial data. In this section we present the tech- 
nical details of the error estimates in the case of ff^(M)-initial data. In this case we do 
not require the numerical scheme to be dispersive, the only ingredient being the Sobolev's 
embedding H^{R) ^ L°°(M). 

Theorem 6.1. Let f{u) = \u\Pu with p G (0,4) and u G C(R, H^{W)) he solution of (fL2]l 
with initial data ip S H^iM). Also assume that is an approximation of order two of the 
laplace operator and u^ is the solution of the following system 

f iu'l + AhU^ = f{u^), t>0, 

(6.3) 

[ ^.'^(O) = Thcp, 

satisfying \\u^\\L^((o,T)xhZ) < C{T,\\ip\\Hi{R)). 

Then for all T > and h > 
(6.4) 

\\u''{t)-Thu{t)\\i2(^hZ) < ^^^^max{r,r2}(||v9||^i(ig) + ||u||^t\o,T://l(K))) '^^P(^II^IIl°°(0,T;H1(IR)))- 

We now give an example where the hypotheses of the above theorem are verified. We 
consider the following NSE: 

, s J iut + d^u = \u\Pu, x G E, t > 0, 

^ ' \ u{0,x) = <pix), X G M, 
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and its numerical approximation 

r iv!l + ^hu^ = W^^vf", t > 0, 

In the case of the continuous problem we have the following conservation laws (see ^ , Corol- 
lary 4.3.4, p. 93): 

lk(*)llL2(M) = II'/^IIl2(R) 

and 

I |n,(t,x)|2dx + — ^- / \uit,x)\P+^dx) =0. 

The same identities apply in the semi-discrete case (it suffices to multiply the equation ()6.6p 

h ■ 

identity) : 
and 



-r(- f \ux(t, x)\'^dx -\ 

api 

by u^, respectively uj*, to sum over the integers and to take the real part of the resulting 



In view of the above identities, the hypotheses of Theorem 16.11 are verified. 



Proof of Theorem \6.1[ Using the variations of constants formula we get 

Thu{t) = Th e^itdx)^ + [ Th exp(i(t - a)dl)f{u{a))da 

Jo 

and ^ 

u''{t) = eMitAh)Thf + / exp(i(t - a)Ah)f{u''{a))da. 

Jo 

Then 

(6.7) err;,(t) := \\u''{t) - Thu{t)\\i2^hZ) 

< II ex.p{itAh)Th(p - Th exp{itdl)ip\\i2(^hz) 

+ ^ II exp(i(t - a) Ah) [fiu'^ia)) - Thf{u{a))) da\\ 

+ / \\eMi{t-(^)Ah)Thf{u{a))-TheM{t-<T)d^x)f{n{a))\\p^hZ)dc7. 
Jo 

Now, applying the error estimates for the linear terms as in (|3.9p with e(l, h) = h^^"^, we get 

(6.8) \\exp{itAh))Th(p -Thexp{itdl)ip\\i2(^hZ) < Th^^^\\ip\\Hi{R)- 

Also, using that f{u) = \u\Pu we have that ||/('u)||hi(m) < ^II^II//i(r) then by (|3.9p we 
get 

(6.9) / II exp(i(t - a)Ah)Thf{u{a))-Th exp(i(t - a)dl)f{u{a))\\i2^hZ)dcT 
Jo 

^ <-^J n II"IIl°°(0,T;H1(R))' 
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Using the /2(/iZ)-stability oiexp{itAh), dE?]), and we obtain 

Now we write - T;,/(u(s)) = /J^(s) + /^(s) where 

I'^is) = f{u\s)) - fiTMs)), I^is) = fiThuis)) - Thfiuis)). 
In the case of we use that / satisfies (|6.ip to get 

\\Ii{s)\\PihZ) < C{\\u''{s)\\^^^^^^ + \\Thu{s)\\P^^j^^^)\\u\s)-Thu{s)\\i2(^f,^) 

< C{\\u^\\lo.^^o,T)xhZ) + ll^lli-{(0,T)xR))ll^''(s) - Thu{s)\\i2(^hZ) 
^ ^II^IIloo(o,T;//1(R))'^^^'*('^)- 

Using the same arguments as in Lemma 15.21 we obtain that 

\\I^is)yikZ)<h\\n{s)C,\^y 

Putting together all the above estimates, for any < t < T we obtain: 
err/j(t) < /i^/^T||(^||j^i(k) + 1^11^00 (o,T; (M)) err/,(o-)dcr 

< h^/^ inax{T,T^}{\\ip\\Hi(R) + lkllit\o,T: hi{m))) + II^IIl-(o,T;//i(r)) ^ err/,(s)ds. 
Applying Gronwall's Lemma we obtain 

(6.10) err;,(t) < /i^/^ max{r, T^} + ||«||it\o,T://i(K))) ^^^P^^H^HL-Co.Ti/fHiR))^' 
The proof is now finished. □ 

6.2. Approximating H''(R), s < 1/2, solutions by smooth ones. Given ip G we 
choose an approximation (p £ H^{^) such that ||v5~'^||h=(r) is small (a similar analysis can be 
done by considering ips G H'^^ with si > 1/2). For ip we consider the following approximation 
of u solution of the NSE (II. 2p with initial data (p: 

idtUh{t) + AhUh = f{uh}, t > 0, 

(6.11) 

Uh{0) = Th'P, 

where the operator Ah is a second order approximation of the Laplace operator. We do not 
require the linear scheme associated to the operator Ah to satisfy uniform dispersive estimates. 

Solving (|6.1ip we obtain an approximation Uh of the solutions u of NSE with initial datum 
(f, which itself is an approximation of the solution u of the NSE with initial datum ip. 

In the following Theorem we give an explicit estimate of the distance between Uh and u. 

Theorem 6.2. Let < s < 1/2, ip E H^{R), and u G C(M; i?'*(M)) he the solution of NSE 
with initial datum ip given by Theorem 15. For any T > there exists a positive constant 
C{T, \\<p\\l^(m.)) such that the following holds 
(6.12) 
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for all h > and 6 > 0. 

In the following we show that the above method of regularizing the initial data (p G (M) 
and then applying the i7^(M) theory for that approximation does not give the same rate of 
convergence /i*/^ obtained in the case of a dispersive method of order two (see (|5.18p ). This 
occurs since for Hv? — <^||L2(ig) to be small, needs to be large and H'ulliooj-Q^ji. 

too. 

To simplify the presentation we will consider the case p = 2. 

Theorem 6.3. Let p = 2, < s < 1/2, (p e H^{R) and u G C(M, H^{M.)) he solution of NSE 
with initial data ip given by Theorem \5.1\ and be the best approximation with (W) -initial 
data as given by (j6.1ip with the conservative approximation = A^. Then for any time T, 
there exists a constant C{\\ip\\H'>{R),T, s) such that 

(6.13) ||TfeU-<||ioo(o,T;/2{/iZ)) <C{\W\\H''{w),T,s)\\ogh\~~-. 

To prove this result we will use in an essential manner the following Lemma. 



Lemma 6.1. Let < s < 1 and h G (0, 1). Then for any ip € H^{M.) the functional Jh^^ 
defined by 

(6.14) Jh,A9) = \\W-9\\h{K) + ^^^v{\\9fmm) 



satisfies: 

(6.15) min Jh,^{g) < C{\\(p\\Hs(R),s)\logh\ 



-s/{l-s) 



Moreover, the above estimate is optimal in the sense that the power of the \ log h\ term cannot 
be improved: for any < e < 1 — s there exists ipe € H^(M.) such that 

(6.16) ^.^.^ min,^^.(^) J.,,.(g) 

^ ^ h^o I log/i|-(^+=)/(i-^-=) 

Remark 6.1. We point out that, to obtain (j6.15p and (j6.16p . we will use in an essential 
manner that s < 1. In fact in the case s = 1 the minimum of Jh over i/^(M) is of order h. 
This can be seen by choosing g = ip and observing that Jhif) = ^exp(||93||j|^i(]j)). This choice 
cannot be done if ^ ^ H'^{E)\H^( 



Proof of Theorem 1 6. 31 Let us choose (p G H'^i^) which approximates (p> in H^{M.). Then by 
Theorem 16.21 we get 

(6.17) \\ThU - ^th|lioo(o,T;«2(hZ)) ^ '^(^' IIv2||l2(r))II<^ - '^lli2(K) + /i exp (2r||(^||2^i(^)) 

< C{T, ||v?||i2(R))4^^^(\/2r(^), 

where Uh is the solution of (j6.1ip with initial data T/j(^. 

For each h fixed, in order to obtain the best approximation u\ of T/,,u, we have to choose 
in the right hand side of the above inequality the function ip* which minimizes the functional 
Jfi ^/2f if (') defined by ()6.14p over H^(U). Using estimate ()6.15p from Lemma [6TT] we obtain 
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the desired result: 



l|T/iU-<||ioo(o,T;/2(hZ)) <C{\\ip\\Hs{M.),T,s) _ TOIU J {V^T (p) 

< C'(ll¥'llH={R),^,'S)|log/ir^ 
where is the solution of (j6.1ip with initial data T^ip* 



□ 



Proof of Lemma \6.1\ The functional J^^^ is convex and its minimizer, g^, is unique. The 
function g^^ satisfies the following equation: 

-^ + 9h + helC£,{\\gh\\]Jll^^){-l^gh + gh) = 



(6.18) 
and so 

I + hexp{\\gh\\Hi^^^)iI - A) 
Thus Ch = \\gh\\H^{R) is the unique solution of 
(6.19) c,= (/-A)V2[/ + /,exp(c2)(/-A)]-' 



9h = (P- 



Step I. A useful auxiliary function. Let us consider the function qh{x) = hx^ exp(x) — c 
for some positive constants /3 and c. We prove that there exist two constants ai(c) and 02(0) 
such that the solution Xh of the equation Qhix) = satisfies 

(6.20) |log/i| - /3 log I log /i| +ai(c) < Xh < \ logh\ - /3 log | log /i| +02(0). 

Let us choose a real number a. Using that h = exp(— | log h\) we get: 

qh{{\ log/i| - /3 log I log/i| + a)) = (| log h\ - /31og | log h\ + af exjp{-^ log | log h\ + a) - c 

log I log /i| ^ Y f \ 

f^+M exp(a)-c 

I log h\ I log h\ J 

-^h^o exp(a) - c. 

Choosing now two constants ai and 02 such that exp(ai) < c < exp(a2) and using that the 
function qf^ is increasing we obtain that, for h small enough, x^, solution of qhix) = 0, satisfies 

Step II. Upper bounds on Ch- Using that ip G H^{R), identity ()6.19p gives us 



Ch 



(j_A)i/2 [l + heMcl){I-A)]-\ 



(/ _ A)(i-)/2 [I + hexpicl)il - A)] ' (/ - A)^/V 



(/ie^')(^-i)/2 r/ie'='(/_A) 



(l-s)/2 



L2( 

-1 



/ + /ie"K/-A) (/-A)^/V 



L2( 



< (/iexp(c;,)) 



,2^^(s-l)/2| 



since, when s G [0, 1], the symbol in the Fourier variable of the operator 



he'^il-A] 



(l-s)/2 r 



I + he^h(i -A) 



is less than one. 
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Then c2(/iexp(c2))i ' < and 
(6.21) {cl)'/^'-'^he< < ll^ll^/.y . 

Applying the result of Step I to /3 = 1/(1 — s) and c = II vll^/'ifR-)'''' obtain that Ch satisfies: 



(6.22) 



cl < I log h\ 



for some constant 02 = a2(||¥'||^iL')**'*)- In particular, when s < 1 



1-s 



log I log /l| + 02, 



/lexp(c^) = exp(c| — | log/i|) < exp 



1-s 



log I log /l| + 02 W 0, 



as /i — )• 0. 



Step III. Estimates on Jhidh)- Using that the minimizer satisfies equation (j6.18p 

and Ch = WghWmiR), we get 

2^^min^^ ^(5) = '^Jhi.gh) = llv' - 5/i||l2{r) + ^ expdl^^H^i^jj)) 
= ihexp{cl)f\\il - A)gh\\l2(^^-j + hexp{cl) 

= (heMcDfUl - A) [/ + heMcl)iI - A)]"' + /iexp(c|) 



l-s/2 



/ + /ie^h(/-A) (/-A)^/2 



-1 



< ihe'^^'Mls^R) + he'^"^ < {he'^ly (||^|||,.(m) + [he^-f'') 

< C(S, ||V9||H=(IR;))(/Iexp(4))^ 

where in the last inequality we used that s < 1 and /iexp(c|) — )• as /i — )• 0. 
Thus, by (l6:22]l we obtain that 



(6.23) 



^J^in^^ ^(5) < c(s, ||99||H-(R))(/iexp(c^))'' < c(s, ||99||/^«(k))| log /i| 1- 



Step IV. A particular function ip. Let us choose e > and (/J^ be defined by means of 
its Fourier transform 



:i+e 



Thus, for any e > 0, <y9e G H^{K). We will prove that, in this case, the solution Ce^h of ()6.19p 
satisfies 

1 



(6.24) 

and 
(6.25) 

for some constant Oi. 



c^h > I log h\ 



1 — S — £ 



log I log h\ + oi 



min Jh,^^{g) > (/iexp(c,\))^+^ > | log /i|-(^+^)/(i-^-^). 
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To prove (j6.24p and (|6.25p we claim that for any 7 G (—1/2,2) and x large enough the 
following holds: 



{x + l + e? ~ x3/2-7- 

Using that C/^^h is solution of (j6.19p and estimate (j6.26p with 7 = 1/2 — s — e and x 
{h exp f^)~^ we obtain 

"^•'^ ^(l + /,exp(c2,)(l + e2))2^^ 

(1 + a--^ 



> 



1 



(/iexp(c2 



and 



/.exp(c2,)(c2,)V(i— )-l>0. 
Applying Step I to the function = /ix^/^-'^^*"^^ exp(x) — 1 we find that 

(6.27) c2 ^ > I log /i| - 3-3737 log I log /i| + ai , 

for some constant oi. 

This concludes the proof of (16.24p . 

We now prove ()6.25p . In view of ()6.18p the minimizer g^^h satisfies 

(6.28) -Vs+ 9e,h + hexp{\\g^^h\\Hi(R)){-'^9e,h + 9e,h) = 
and 

(6.29) ge,h= [/ + /iexp(c2 j(/- A)]"Ve, 

where Ce^h = \\ge,h\\HHR)- 
Thus 

'^Jh,<fi,{ge,h) = \\(Pe - gsAWh^tL) + ^^^Pi\\9e,h\\m{M.)) 

= {hexp{clf^)f\\{I - A)gs^h\\l2(^^-j +hexp{clf,) 

= ihexp{clh)f\\il - A) [l + /iexp(c2 J(/- A)]"Ve|li2(K) + /iexp(c2 J. 
Writing the last term in Fourier variable we get 



2 mm Jh,^^{g) = {hexp{c^f^)) / + hexp{c,^h 

-^^ (l + /iexpfc2,)(l + £2)) 



^ +/iexp(c2j(l + e2) 

^-^^ + hexp{cl„). 

{hexp{cl,))-^ + l + e) 
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The same arguments as in Step II give us that /iexp(c^^) — )■ as /i — t- 0. Then for small 
enough h, Xh defined by = (/iexp(c^ ^))~^ is sufficiently large to apply inequality ()6.26p 
with 7 = 3/2 — s — e. We get 

2 "^Y}^. '^'^'^^^S) > iheMclh)y^' + hexpiclf,) > (/i exp(c2^))^+^ 
Using now ()6.27p we obtain 



min J/, (5) > I log /i| 

which proves (j6.16p . 

To finish the proof it remains to prove (j6.26p . For |x| — >• c«, using changes of variables we 
get: 



00 27 



^=xi/2£ j;3/2-7 y^>^._i/2 1 + e 

+ 0(-^)>^7^(l + 



> 



J.O/Z— 7 ■ ' ~ 2;3/2— 7 ^ ^1/2+7' 

1 



which proves ()6.26p . □ 
Proof of Theorem \6.^ Using the stability result (|5.ip for the NSE we obtain 

\\u - 'u||l°°(0,T;L2(R)) < CiT,P, ||V^||l2(]j), ||<^||l2(r))||(^ - 
< C{T,p, ||<^||z,2(R))||(^ - 1^||l2(k). 

Now using the classical results for smooth initial data presented in Section 16.11 by (j6.10p 
we get 

W^hU - Uh\\L^{o,T;i^{hi.)) < C'/i^/^exp(r||n||^^^Q^y.^i^jg^p. 

Thus 

\\ThU - ■U/i||l°o(0,T;/2(/iZ)) < W^hU - T/iU||ioo(o^T; /2(/iZ)) + llTftU - Uh\\L°°(0,T;P{hZ)) 

< II'" - '"IIl°°{0,T;L2(K)) + ||T/itt - n/i||L°°(0,T;«2(/iZ)) 

< C(r,p, ||(^||i2(K))||v3 - (^||l2(k) + /i^/^ exp(r||i2||^^(Q^y.^i(]g)P. 
This yields (I6T2D. □ 
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7. Technical Lemmas 

In this section we prove some technical results that have been used along the paper. The 
main aim of this section is to obtain estimates on the difference f(Thu) — Thf{u) in auxihary 
norms Li{I,F{hZ)). 

In the case of smooth enough functions u, the pointwise projection operator 

(7.1) {Bhu){jh)=u{jh) 

makes sense. More precisely it is well defined in H^{M), s > 1/2. In these cases the use of 
the operator E/i has the advantage of commuting with the nonlinearity f(Ehu) = 'Ehf{u). 
The key ingredient is the following Theorem. 

Theorem 7.1. (Marcinkiewicz multiplier theorem [13j, Th. 5.2.2, p. 356) Let m : M — t- M 6e a 
hounded function which is in every dyadic set (2^ ,2^^^)U {—2^'^^ , —2^) for j G Z. Assume 
that the derivative m' of m satisfies 



(7.2) sup 



_2i 

\m'{i)\di+ \m'{i)\di 



< A <oo. 



2J+1 J 23 

Then there exists a positive constant C such that for all 1 < q < oo the following holds: 
\\{fm)'^\\Li(R) < Cmax{q, {q - 1)~^}^{A + \\m\\Loom)\\f\\Li(R)- 



Remark 7.1. Using a change of variables in the Fourier space the above dyadic intervals can 
be replaced by any other one of the form {c^,c^^^) U {—c^'^^,—c^), j G Z and c > 1. In the 
following applications, the constant c will be chosen to be c = it. 

For any function u G L^(M) we define the new function Hh by truncating the Fourier 
transform as follows: 

(7.3) UhiO = U{0 l(-7r/h,7r/h)(0- 

For h = 1, Theorem 17.11 applied with m(^) = 1(_^ ,^) which is in every dyadic interval, 
shows that for any 1 < q < oo, the L''(M)-norm of ui can be controlled by the one of u: 



(7.4) \\ui\\LHR)<Ciq)\\u\ 



A scaling argument shows us that the above inequality also holds for all /i > with a constant 
C{q) independent of h. 

Using Theorem 17.11 we can refine this estimate as follows: 

Lemma 7.1. For any s > and q G (l,oo) the following hold, 
a) There exists a positive constant c{s,q) such that 



(7.5) \\u-iih\\Li(R) <cis,q)h^\\u\ 



holds for all u £ Ty^''?(M) andh>0. 

b) Assuming s G [0,1], there exists a positive constant c{s,q) such that 
(7-6) h\\uh\\^,i,g^^) <c{s,q)h'\\u\\y^,s:,^^^ 

holds for all u G VF'''«(M) andh>0. 
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Proof of Lemma 7.1 We divide the proof in two steps corresponding to ()7.5p and (|7.6p . 
Step I. Proof of ()7.5p . Let us consider the fohowing operator 

Mhu ■=u-Uh = {l{\^\>T,/h)uY ■ 
A change of variables gives us that 

(M^n)(x) = Afi(n(/i-))(^). 
Using this property the following identities hold: 

||Mfeu|U,(M) = /ii/'?||Mi(n(/i-))llL«(R) 

and 

ll«(^-)lliy=..(R) = ll|Vr[^x(/i-)]IU.(R) = /i1|(|Vrn)(/i-)||L,(M) = h'h^''^\\\V\'u\\L,i^^y 

Thus, it is sufficient to consider the case h = 1 and to prove that 

(7.7) <c(s,g)|||V|'*^;||i,(K) 

holds for aU v £ W^'HR). 
With the notation 

m,(0 := l^r^l{|C|>.}(0, 
estimate (|7.7p holds if ms{S,) satisfies the hypothesis of Theorem 17.11 Using that ms(^) G 

<(e)i<^i{i5i>^}(0,eGK, 



L°°(R) and that 

\m' (^M < 



by Theorem 17.11 we obtain (j7.7p . 

Step II. Proof of (|7.6p . A similar argument as in the previous case reduces estimate 
to the case h = 1: 



ii(s(oi(_.,.)iei)^iiL.(M) < c{s,q)\\{um\r\\LHm- 

Denoting v = (u{^)\^\^)'^ , it remains to prove that 

(7.8) ll(?^(Ol(-.,.)ler-')''llL.(M) < c{s,q)\\v\\L.iR). 

In other words, it is sufficient to apply Theorem 17.11 to the multiplier ms{(,) given by 

m,(e) = iei'-^i(-.,.)(0- 

Using that ms{0 ^ i°°(K) satisfies 

im',(e)i <c(5)ier^i(-.,.)(6,^eM\{o}, 

we fit in the hypothesis of Theorem 17.11 and then (j7.8p holds. This finishes the proof. □ 

In the following we obtain error estimates for the difference between the two interpolators 
T/j and when applied to functions u and f{u), where and E/j are defined by (j3.2p and 
(|7.ip respectively. 

Lemma 7.2. Let s > 1/2 and q £ (l,oo). Then there exists a positive constant c{s,q) such 
that 

(7.9) \\ThU-Ehu\\ig(_hZ) < c{s,q)h''\\u\\^s,g(^^) 
holds for all u e W'^iR) and h>0. 
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Remark 7.2. This lemma generalizes Theorem 10.1.3 of [30j, p. 205, which addresses the 
case q = 2, s > 1/2. In this case using Plancherel's identity in the discrete setting it is easy 
to obtain 



(7.10) 



\ThU-Ehu\\i2(^h^) < c{s)h''\\u\\^s 



Remark 7.3. Using the above results, we will be able to obtain estimates of the difference 
T^hfiu) — f(Thu), f{u) = \u\Pu, p> 0, given by Lemma [57B . 



Proof of Lemma 7.2. Estimate (17.10p provides the desired estimate VF*'^(]R) — ?■ P{h1i) in the 
case q = 2. We will also prove the estimate W^''' — )• l'^{h'Z) in the case s > 1. Using these two 
estimates the general case will be a consequence of an interpolation argument. 
Case 1: s > 1, q e (l,oo). We claim that 

(7.11) \\ThU-Bhu\\ip^hZ) <c{p,s)h'\\\V\'u\\LP(R). 
By rescaling all the above quantities we can assume h = 1. 

We have the following: 

(Tiu-Ein)(j) = / e'^^u{i)= I e'^^^u{i + 2^1). 

J\(,\>1T J-TT 

Denoting by v the function whose Fourier transform is given by 

(7.12) v{0 = li-^,n)Y.u{^ + 27:l), 

we get 

(Ti^x-Ei^x)(i) = r e'^^v{i)di. 

J — TT 

Classical results on band-limited functions (see Plancherel-Polya [27j) give us that 

\\Tiu -Biu\\ip(^z) < II^IIlp(R)7 

provided that the right hand side term of the above inequality makes sense. It is then sufficient 
to prove that the function v defined by ()7.12p satisfies: 

(7.13) IbllLp(R) < c(p,s)|||V|''n||ip(i;). 
Using the properties of the Fourier transform we get: 

v{x) = J^e2^^'^(l((2z_i)^,(a+i)^)n)'^. 

It is sufficient to prove that 



e^''''^(l((2/_l),r,(2/+l)7r)^^) 



< \\M"u\ 



\LP( 



LP( 



or equivalently 



-'-((2/ — l)7r,(2/+l)7r)^J 

1^0 



< u 



LP( 
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We claim that for any / 7^ 0: 

(7.14) ||(ier'l{(2^-iK(2mwn)^| 



< '^'\\ II 

LP - -mTimiLn 



Thus, summing all the above inequalities for I ^ we obtain the desired estimate. 
A translation in (|7.14p reduces its proof to show that m^,;, defined by 

"^s,«(o = ie-2/^r^i(-.,.)(o, 

verify the hypothesis of Proposition 17.11 Observe that 



and 



CIS 



|m',/OI < e G K \ {0}, / / 0. 



Applying Proposition 17.11 to each multiplier nis^i we get ()7.14p and the proof of this case is 
finished. 

Case 2: s > 1/2, q G (l,oo). We set Uh = Th — Ei^. Using the estimates of the previous 
case we deduce that the operator Uh satisfies: 

Uh : W''^'^^{R) I'^'ihZ), si > 1, 1< < 00, 

and by fTTOD : 



Uh-.W 



52,2/ 



l^{hZ), S2 > 1/2. 



Then for any 9 G (0, 1), 

Uh : [W''''i^iR),W'''\R)\e] ^ [I'^^h^) , l\hZ)\e] 
with a norm that satisfies: 

\\Uh\\[wnm(R),W''2>2(R)]^g^-[lii(hZ),P{hZ)]ig] ^ \\^f^\\w^i-n (R)-i'ii (hZ)^^^f'^^\V^2,2(M,)-p(hZ)- 
Classical results on interpolation theory ([2] ,Th. 6.4.5, p. 153) give us that 

and 

l[e] 

where s and q are given by 

s = fe{si, S2) = 516! + S2(l - 0), 



[/^i(/lZ),/2(/iZ)]r0i = l1{hZ) 



(7.15) 



1 / N ^ 1- 

- = 9e{qi) = — + 

^ q gi 2 



Using that the ranks of functions fg and gg satisfy 

■1 + 



Im{fg 



00 , Im{gg) 



1-6 1 + 



2 ' 2 /' 
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we obtain that for any s > 1/2 and < g < 1 we can find Si > 1, S2 > 1/2, qi > 1 and 
9 e (0, 1) such that (f77T5]) holds and 

(7.16) IIA|lvy..,(M)_z,(,z) < h^'%^-('-'^ < h^. 

The proof is now finished. □ 
Proof of Lenima \5.2l We first recall that the following inequality holds for all u,v £ L^'"^^(R): 

(7.17) ll/(^i) - /(i')IIl{p+2)'(k) < C{p){\\u\\%+^ + \\v\\Ip+2)\\u - v\\Lp+2^f,). 

We set Uh defined by Uh{£.) = u{0 l{-7r/fe, 7r/fe)(0- The difference Thf{u) - f{Thu) in (|5.25p 
satisfies: 

l|Th/(u)-/(T/ju)||;(j,+2)'(;j2) < \\Thf{u)-Thf{uh)\\i(p+2y(^hZ) + \\^hf{uh)-f{Thu)\\i^p+2y(^f^j^y 
Using (j7.17p . (|7.4p and Lemma |7. 11 the first term in the right hand side satisfies 
\\Thfiu) -Thf{uh)\\i(p+2y^l,j^) < c{p)\\f{u) - f{uh)\\L(p+2y^^) 

< c{p){\\u\\Ip+2(^^^ + \\Uh\\lp+2(^si)^\\u - %||lp+2(r) 

< c(p)/l*||tl||^p+2(R)||^^|lvi/a,p+2(R) < C(p)/l* II n||^s_p+2(lR)- 

For the second term, using that on the grid hZ, T^u = E^Uh, by Lemma ()7.2p we get: 

(7.18) ||T;,/(n;,)-/(Tftn)||,(,+2/(^^) = \\Thf{uh) - f{^hUh)\\np+2y ^^.^^ 

= \\'^hf{uh)-'E,hfiuh)\\i(p+2y(^hi.) 

< h\\fiuh)\\wi.(p+2y(K) ^ H^h'^xUhW i^(p+2y (^y 

Using that s G [0, 1] we apply Young's inequality and ()7.6p to obtain: 

(7.19) ||<5.n.||,(,,2)/(...)(^) = ( |j^.r(^+^)/(-+^)|a.^.|(^-^^)/(--^^))^'^'^^^'^'^ 

= ll^/i|ll,P+2(R) l|f^a:''Uh||LP+2(R) < ll^ll /^p+2 H'Uh || ^yi,p+2 

Thus by fTTSl) and ([7J9]1 we obtain 

\\T^hf{Uh) - /(Tfe-u)||;(p+2)'(;^^) < /i''||'U||^S^p+2(IR) 

which finishes the proof. □ 
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